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Abstract: We compute exact values respectively bounds of "distances" - in the sense of (transforms 
of) power divergences and relative entropy - between two discrete-time Galton- Watson branching pro- 
cesses with immigration GWI for which the offspring as well as the immigration is arbitrarily Poisson- 
distributed (leading to arbitrary type of criticality). fmplications for asymptotic distinguishability be- 
haviour in terms of contiguity and entire separation of the involved GWI are given, too. Furthermore, 
C_^ , we determine the corresponding limit quantities for the context in which the two GWI converge to 

Cm ' Feller-type branching diffusion processes, as the time-lags between observations tend to zero. Some ap- 

plications to (static random environment like) Bayesian decision making and Neyman-Pearson testing 
are presented as well. 
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1. Introduction 



^- ' It is well known that "distances" in form of (relative-entropy covering) power divergences between finite 

measures are important for probability theory and statistics as well as their applications to various different 
research fields such as physics, information theory, econometrics, biology, speech and image recognition, 
transportation of (sorts of) "mass", etc. For probability measures P-^, P4 on a measurable space (Q,A) and 
parameter A G R these power divergences - also known as Cressie-Read measures respectively generalized 
cross-entropy family - are defined as (see e.g. Liese and Vajda [47], [48]) 



f HPaWPu), if A = l, 

Ix(Pa\\Ph) ■= \ x^T)(Hx(PA\\Pn)-l), if AeR\{0,l}, 

I(Ph\\Pa), if A = 0, (1) 

£■ ™ e I{Pa\\Ph) ■■= I PAtog— d/i + 00 -P a (ph = 0) (2) 

J{pn>0} PH 

is the relative entropy (Kullback-Leibler information divergence) and 

H x (P A \\Pu) := f P A Pn X dfi (3) 

Jn 

is the Hellinger integral of order A € R\{0, 1}; for this, we assume as usual without loss of generality that 
the probability measures P-u, Pa are dominated by some a— finite measure /i, with densities 

dP A , dPu 

PA = -5 — and p-u = -5 — 

defined on CI (the zeros of p-u, p A are handled in (2), (3) with the usual conventions). Apart from the relative 
entropy, other prominent examples of power divergences are the squared Hellinger distance \ Jx/2 {Pa\\Ph) 
and Pearson's \ 2 ~ divergence 2 12 (Pa\\Ph)- Extensive studies about basic and advanced general facts on 
power divergences, Hellinger integrals and the related Renyi divergences of order A £ K\{0, 1} 
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Rx (Pa\\Ph) ■= johr) logHx ^H^ ' with log ° = ~°°' 

can be found e.g. in Liese and Vajda [47], [48], Jacod and Shiryaev [29]. For instance, the integrals in (2) and 
(3) do not depend on the choice of fj,. As far as finiteness is concerned, for A G]0, 1[ one gets the rudimentary 
bounds i 

< Ix(Pa\\Ph) < A(A1) , (4) 

where the lower bound is achieved if and only if P4 = P<n , and the upper bound is achieved if and only if 
Pa-LP-h (singularity). For A <^]0, 1[, the power divergences I\ (Pa\\Ph) and Hellinger integrals H\ (Pa\\Ph) 
might be infinite, depending on the particular setup. For the sake of brevity, we only deal here with the case 
A G [0, 1]; the case A ^ [0, 1] will appear elsewhere. 

Apart from the extensive literature on the relative-entropy cases A(l — A) = 0, for A(l — A) 7^ the evaluation 
of power divergences I\ - respectively their straightforward transforms such as Hellinger integrals H\ and 
Renyi divergences R\ - have been investigated for various different contexts of (probability distributions 
associated with) stochastic processes, such as processes with independent increments (see e.g. Newman [57], 
Liese [44], Memin and Shiryaev [55], Jacod and Shiryaev [29], Liese and Vajda [47], Linkov and Shevlyakov 
[53]), Poisson point processes (see e.g. Liese [45], Jacod and Shiryaev [29], Liese and Vajda [47]), diffusion 
prcoesses respectively solutions of stochastic differential equations with continuous paths (see e.g. Kabanov 
et al. [33], Liese [46], Jacod and Shiryaev [29], Liese and Vajda [47], Vajda [69], Stummer [64], [65], Stummer 
and Vajda [67]): further related literature can be found e.g. in references of the abovementioned papers and 
books. 

Another important class of time-dynamic models is given by discrete-time branching processes, in particluar 
Galton- Watson processes without immigration GW respectively with immigration GWI, which have numerous 
applications in biotechnology, population genetics, internet traffic research, clinical trials, asset price modelling 
and derivative pricing. (Transforms of) Power divergences have been used for supercritical Galton- Watson 
processes without immigration SUPGW for instance as follows: Feigin and Passy [15] study the problem 
to find an offspring distribution which is closest (in terms of relative entropy type distance) to the original 
offspring distribution and under which ultimate extinction is certain. Furthermore, Mordecki [56] gives an 
equivalent characterization for the stable convergence of the corresponding log-likelihood process to a mixed 
Gaussian limit, in terms of conditions on Hellinger integrals of the involved offspring laws. Moreover, Sriram 
and Vidyashankar [62] study the properties of offspring-distribution-parameters which minimize the squared 
Hellinger distance \ I\ii between the model offspring distribution and the corresponding non-parametric 
maximum likelihood estimator of Guttorp [19]. For the setup of GWI with Poisson offspring and nonstochastic 
immigration of constant value 1, Linkov and Lunyova [52] investigate the asymptotics of Hellinger integrals 
in order to deduce large deviation assertions in hypotheses testing problems. 

In contrast to the abovementioned contexts, this paper pursues the following main goals: 

(MG1) for any time horizon and any criticality scenario, to compute (non-rudimentary) lower and upper bounds 
- and sometimes even exact values - of the Hellinger integrals H\ {PaWPh) an d power divergences 
I\ (Pa\\Ph) (A G [0, 1]) of two Galton- Watson branching processes P_a, P-n with Poisson(/3_4) respec- 
tively Poisson(/?^) distributed offspring as well as Poisson(a^) respectively Poisson(a-H) distributed 
immigration. As a side effect, we also aim for corresponding asymptotic distinguishability results in 
terms of contiguity and entire separation. 

(MG2) to compute the corresponding limit quantities for the context in which (a proper rescalation of) the 
two Galton- Watson processes with immigration converge to Feller-type branching diffusion processes, 
as the time-lags between the generation-size observations tend to zero. 

(MG3) as an exemplary field of application, to indicate how to use the results of (MG1) for Bayesian decision 
making and Neyman-Pearson testing based on the sample path observations of the GWI-generation 
sizes, when the hypothesis law is given by Pu and the alternative law by P4; in a certain sense, this 
can also be interpreted in terms of a rudimentary static random environment. 

Because of the involved Poisson distributions, these goals (which are potentially reasonable also for other 
types of offspring resp. immigration distributions) can be tackled with a high degree of tractability, which is 



Niels B. Kammerer and Wolfgang Stummer/ Distances of branching processes 3 

worked out in detail with the following structure: we first deal with the non-relative-entropy case A(l — A) 7^ 0. 
Section 2 contains the first basic result concerning Goal (MG1), which is then deepened in Section 3 in order 
to obtain - parameter constellation dependent - recursively computable exact values respectively recursively 
computable lower and upper bounds of H\ [P^\\P%). Additionally we construct related closed-form bounds 
in Section 4, which will also be used to achieve (the Hellinger-integral part of) Goal (MG2) in Section 5. 
The power divergences I\ {Pa\\Ph) & re treated in Section 6, complemented with the relative-entropy cases 
A(l - A) = of the Goals (MG1), (MG2). The subsequent Section 7 is concerned with Goal (MG3), whereas 
the Appendix contains main proofs and auxiliary lemmas. 

2. Process setup and first basic result 

Let X n denote the nth generation size of a discrete-time Galton- Watson process with immigration GWI. We 
use the recursive description 

X»-i 

A := uo a G N; X n = £ F„_ 1>fc + Y n , n e N, (5) 

fc=i 
where Y n —ik is the number of offspring of the fcth object (e.g. organism, person) within the (n — l)th 
generation, and Y n denotes the number of immigrating objects in the nth generation. Notice that we employ an 
arbitrary deterministic initial generation size Xq. We always assume that under the law Ph (e.g. a hypothesis), 

• the collection Y := {Y n —i,k: n € N, fc € N} consists of independent and identically distributed (i.i.d.) 
random variables which are Poisson distributed with parameter Ph > 0, 

• the collection Y := \ Y n , ngN consists of i.i.d. random variables which are Poisson distributed with 
parameter a-n > (where au = stands for the degenerate case of having no immigration), 

• Y and Y are independent. 

In contrast, under the law P4 (e.g. an alternative) the same is supposed to hold with parameters 0a > 
(instead of Ph > 0) and a a > (instead of a% > 0). Furthermore, let (J 7 n)ne^ be the corresponding 
canonical filtration generated by A := (A"„)„ 6 n- 

Basic and advanced facts on GWI (introduced by Heathcote [21]) can be found e.g. in the monographs of 
Athreya and Ney [2], Jagers [30], Asmussen and Hering [3], Haccou [20]; see also e.g. Heyde and Seneta [25], 
[25], Basawa and Rao [4], Basawa and Scott [6], Sankaranarayanan [59], Wei and Winnicki [71], Winnicki 
[72], Guttorp [19] as well as Yanev [73] (and also the references therein all those) for adjacent fundamental 
statistical issues including the involved technical respectively conceptual challenges. 

For the sake of brevity, wherever we introduce or discuss corresponding quantities simultaneously for both the 
hypothesis H and the alternative A, we will use the subscript • as a synonym for either the symbol % or A. For 
illustration, recall the well-known fact that the corresponding conditional probabilities P. (A„ = • |A„_i = fc) 
are again Poisson-distributed, with parameter 0, ■ k + a,. In oder to achieve a transparently representable 
structure of our results, we subsume the involved parameters as follows: let Psp be the set of all constella- 
tions (0a> Pv.1 a Ai ol<h) of real- valued parameters Pa > 0, 0h > 0, cva > 0, a-u > 0, such that Pa 7^ Ph or 
a A 7^ a H ( or both). Furthermore, we write Pni for the set of all (Pai Ph-, a A-, a u) of real- valued parameters 
Pa > 0, Pu > 0, cia = a u = 0, such that Pa ^ Ph] this corresponds to the important special case of 
having no immigration. The resulting disjoint union will be denoted by V — Psp U Pni- A typical situation 
for applications in our mind is that one particular constellation (Pa> Put &Ai ecu) €E V (e.g. obtained from 
theoretical or previous statistical investigations) is fixed, whereas - in contrast - the parameter A s]0, 1[ for 
the Hellinger integral or the power divergence might be chosen freely, e.g. depending on which "probability 
distance" one decides to choose for further analysis. At this point, let us emphasize that in general we will 
not make assumptions of the form /3. = 1, i.e. upon the type of criticality. 

To start with our investigations, we define the extinction time r := min{Z G N : X m = for all integers 
m > 1} if this minimum exists, and r :— 00 else. Correspondingly, let B := {r < 00} be the extinction set. 
It is well known that in the case Vm one gets P,(B) — 1 if < 0, < 1 and P,(B) G]0, 1[ if 0, > 1. In 



Niels B. Kammerer and Wolfgang Stummer/ Distances of branching processes 



contrast, for Psp there always holds P,{B) = 0. Furthermore, for Vsp the two laws P% and P4 are equivalent, 
whereas for "Pni the two restrictions Ph\b an< ^ -^Uls are equivalent (see e.g. Lemma 1.1.3 of Guttorp [19]); 
with a slight abuse of noation we shall henceforth omit \ B . Consistently, for fixed time n G No we introduce 
Pa,u '■= Pa\jt an d P-u,n ~ Ph\jt as we ll as the corresponding Radon-Nikodym-derivative 

dPA,n 



dP n ,; 



(6) 



Clearly, Zq = 1. By using the "rate functions" f,(x) = (3,x + a, (x G [0, oo[), a version of (6) can be easily 
determined by calculating for each u> = (wq, • ■-, w n ) G &n '■= Nq 



n 

Z n (u) = \\ Z n . k {w) with Z n . k {uj) := exp | - (/^(wfe_i) - /^(wfc-i)) j 



fe=i 



where for the last term we use the convention 



JU(w fe - 



/«(Wfc-l) 



'O^ 1 



CJ G ilr 



?$(w) := exp { - (A/ACwfc-O + (1 - A)/ W (w fc _i))} 



(/^(Wfc-l)) (/•H(wfe-i)) 1 



1 for all x G No- Furthermore, we define for each 

wfc 

(7) 

with the convention M- = 1 for the last term. Accordingly with the choice /1 = Pu,n one obtains from (3) 
the Hellinger integral H\ (Pa,o\\Pv.,o) = 1, as well as for all (j3a,P~h,&a, a-u. A) € ~Px]0, 1[ 

i?A (Pa,iII^,i) = cx p{ (/-aM) A (/«M) (1 ~ A) - (A/a(w ) + (1 - A)/«(wo))} 
and for all n G N\{1} 



(8) 



^(PAnll^.n) - EPn,n[(Z n ) X ] = £ • • • £ JJ Zg> (c 



e- e no 

CJl— CJn-1— Ofc — 1 

00 00 n — 1 



^1=0 w„=0fc=l 



Wf w ) . e -(A/^(w„_i)+(l-A)/K(u; B _i)) 



E 



[(/yt(^„-l)) A (/-H^-i)) 1 " 



e- e n*-« 



(/^(w„_i)) A (/«( W „-i)) 1 



>=0 



-(A/x(w«-i)+(1-A)/m(w„_i)) 



(9) 



wi- u n -i— fc— 1 



From (9), one can see that a crucial role for the exact calculation (respectively the derivation of bounds) of 
the Hellinger integral is played by the functions defined for x G [0,oo[ 

<j>\(x) := ifi\(x) - f\{x) , with (10) 

<P\(x) := (f A (x)) X (f n (x)) 1 - X and (11) 

fx(x) := Xf A (x) + (l-X)Mx) = a x + p x x, (12) 

where we have used the X-weighted-averages )3\ = A • /3„4 + (1 — A) ■ /3-n and «a = A • a A + (1 — A) • a%. 

According to Lemma A.l in Appendix A.l, it follows for A G]0, 1[ that <f)\(x) < for all x G [0, oo[, and that 

4>{x) = iff fA( x ) = fu(%)- This is consistent with the corresponding generally valid upper bound 

Hx(PA,n\\PH,n)<l. (13) 

As a first indication for our proposed method, let us start by illuminating the simplest case A G]0, 1[ and 
7 := ch-uPa ~~ a A^H — 0. This means that (/3ajPh-> OL A^ a 'K) G Pni U Psp.Ij where Psp.i is the set of all 
(componentwise) strictly positive (Pa,Ph,&a,&h) with f3 A ^ (3-u, a a 7^ cx-h and 4p- = ^ 7^ 1. In this 
situation, aH the three functions (10) to (12) are linear. Indeed, 

fx(x) = pf + qfx 
with pf := a\ a]^ and qf := P\P]7 (where the index E stands for exact linearity). Clearly, qf > on 
Pni U Psp,i, as weu a s pf > on Psp.i respectively pf = on Pni. Furthermore, 

A (z) = rf + sfz 

with rf := pf - oa = a* a^ A - (Aoa + (1 - A)o w ) a nd sf := qf~Px = P A P X U X ~ (X(3 A + (1 - X)f3 n ). 
Due to Lemma A.l one knows sf < on 'Pni UPsp,i , as well as rf < on Psp,! respectively rf = on Pni- 
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As it will be seen later on, such kind of linearity properties are useful for the recursive handling of the Hellinger 
integrals. However, only on the parameter set Pni U "Psp,i the functions ip x and <f> x are linear. Hence, in the 
general case (/3a, /3-h, a a, a-u, A) G ?x]0, 1[ we aim for linear lower and upper bounds 

<P\(x) -=P\ +Q\X < <p\(x) < <p%(x) ~P\+q\X , (14) 

x G [0, oo [ (ultimately, x G No), which lead to 

^(x):=rj: + sj:x-.= (pi-a x ) + (q^-p x )x< ( l }x (x)<^(x):=r^ + s^x:=(p^-a x ) + (q^-p x )x, (15) 

x G [0,oo [ (ultimately, x G No). Of course, the involved slopes and intercepts should satisfy reasonable 
restrictions. For instance, because of the nonnegativity of ip x we require p^ > p^ > 0, q^ > q^ > (leading 
to the nonnegativity of p\, v?a )• Furthermore, (9) and (13) suggest that p\ < a X) q^ < f3 x which leads to 
the nonpositivity of c\>\ . Moreover, it is assumed that 

at least one of the two inequalities p\ < a x , q x < [3 X holds, (16) 

and hence 4> x (x) < for some (but not necessarily all) x G [0, oo[. Notice that in (16) we do not demand 
the validity of both inequalities, which might lead to the effect that the constructed Hellinger integral upper 
bounds have to be cut off at 1 for some (but not all) observation horizons n G N; see (21) below. For the 
formulation of our first assertions on Hellinger integrals, we make use of the following notation: 

Definition 2.1. For all (/3a, /3u, a Ai a 'H, A) G Px]0, 1[ and all p G [0,oo[, q G [0,oo[, let us define the 
sequences ( a„ J and ( bn ) recursively by 

\ / nPNn V /nGNn 



a { o q) :=0 ■ a$ := e<-> ■ q - f3 X , n G N, (17) 

b^' q) :=0 ; b%'d := e a --i ■ p - a x , n G N. (18) 



Notice the interrelation a[ 9x ' = s x and 6^ ' qx — r x for A G {E, L, U}. Clearly, for q g]0, oo[, p G [0, oo[, 

one has the linear interrelation , ., v < \ v 

&<M) = £ ai 9) + V - (3 X - a x , neN. (19) 

Accordingly, we obtain fundamental Hellinger integral evaluations: 

Theorem 2.2. (a) For all (/3_4, f3-n,aA,ccHj A) G ( - P./v/U'Psp i i)x]0, 1[, all initial population sizes u>o € N and 
all observation horizons n G N one can recursively compute the exact value 

n 



fc=i 



where %^ can be equivalently replaced by |p-. Recall that qf := P\P H 

(b)Forail(PA,0H,aA,<XH,\) G (^5p\^SP,j)x]0, 1[, all coefficients p^ g [0,oo[, o^ G [0,oo[,p£ G [0,min{p^,a A }], 
g^ G [0, minlg^, /3a}], smc/i i/iai (14) ZioZds /or a// X G No as lueZZ as (16), a// initial population sizes 
wo G N and all observation horizons n G N one geis i/ie recursive (i.e. recursively computable) bounds 
Bi n < H X {PA,n\\ P n,n) < B^ n , where 



B X,n 



:= ex P {a { / x) u Ja + J2b { k pU " ) } and B% n := min exp {a^ co + £ bf'^} , 1 \ . (21) 

Remark 2.3. From the proof below one can see that both parts of Theorem 2.2 remain true for the cases 
A ^ [0,1]. For the (to our context) incompatible setup of GWI with Poisson offspring but nonstochastic 
immigration of constant value 1, the exact values of the corresponding Hellinger integrals (i.e. an "analogue" 
of part (a)) was established in Linkov and Lunyova [52]. 

Proof: 

We first prove the upper bound B\ n . Let us fix ((3a, Ph, a A, a Hi A), p\ , q x , ujq G N as described in part (b). 
From (8), (10), (11), (12) and (14) one gets immediately B x 7 lJ and with the help of (9) for all observation 
horizons n G N\{1} (with the obvious shortcut for n — 2) 
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oo oo n—X 

H x (P A , n \\P n ,n) = E'" E l[ Z^) ■ cxp {^(0J n ^) - f x (0J n ^)} 

UJ\— LU n -i— fe=l 

oo oo n— 1 

< E-- E II z Sh- cx p{(pa-«a) + (^-/3a)^- 1 } 

LJl — OJ„_i— fc— 1 



Ji— cj^-i— fe=l 

CXD 71 — 1 



OO OO 71 — 1 

e- e n^ A >)-^K^ ) +«^ ) «-i} 

^1=0 W n _i=0 fe=l 



oo oo n — 2 

= e XP {^ M } E--- E n z 2H' ex p{ exp { a ^V AK - 2)_/AK - 2) } 

Wl- LO n -2— fc=l 

OO CO 71 — 2 

<exp{^} E- E IK^") 

UJ\— CJ n _2=0 fc=l 

• exp { (exp {a^ } p u x - a x ) + (cxp [a[ q " ] } <# - &) • w„_ 2 } 

OO CO ra — 2 

<exp{^} E- E n*)-«P{^+#^} 

LJi— HJ T1 _2— fc=l 

n 

< ... < expja^^o + E 6 ^^}- ( 22 ) 

fc=i 

Notice that for the strictness of the above inequalities we have used the fact that <f>\ (x) < A (x) for some 
(in fact, all but at most two) x € No (cf. (p-xiv) below). Since for some admissible choices of P\,<1\ and 
some n £ N the last term in (22) can become larger than 1, one needs to take into account the cutoff-point 
1 arising from (13). Notice that without assumption (16), the last term in (22) would always be larger than 
1 (and thus useless). The lower bound B\ n of part (b), as well as the exact value of part (a) follow from (9) 
in an analoguous manner by employing p x , q x and pf , qf respectively. Furthermore, we use the fact that for 

(/3 A ,p H ,a A ,a n ,\) e CPniU:P S p,i)x]0, 1[ one gets from (19) the relation b^' 9 ^ = f^a { ^ ] . For the sake 
of brevity, the corresponding straightforward details are omitted here. Although we take the minimum of the 
upper bound derived in (22) and 1, the inequality B Xn < B Xn is nevertheless valid: the reason is that for 
constituting a lower bound, the parameters p A ,<Z A must fulfil either the conditions [p x < and q x < 0] or 
\p{ < and q% < 0] (or both). □ 



3. Detailed analyses 

For part (b) in Theorem 2.2, we have assumed the existence of reasonable linear lower and upper bounds 
of tpx and cj>\. In the following, we shall carry out a more detailed analysis addressing questions upon the 
non-uniqueness (and thus, flexibility) of the coefficients p x , q x , p^, q x in (14), their "optimal respectively 
reasonable choices", as well as the corresponding behaviour of the Hellinger integrals H\(P Atn \\Pn >n ) as the 
observation horizon n increases and finally converges to oo. Of course, the answers to these questions will 
depend on the (e.g. fixed) value of (f3 A , /3%, a A , an) and the (e.g. selectable) value of A. 

Before starting a closer inspection, notice by induction the general fact that for {j3 A , /3-n, a A , a-u, A) ePx]0,l[ 
and q g]0, oof the principal behaviour of the sequence ( a„ J is strongly governed by its first element: 

V /nGN 

(p-i) a { n q) = 0, if a[ q) = q - /3 A = (i.e. q = X ), 



(p-ii) ( a„ ) is strictly negative and strictly decreasing, if af < 0, 

(p-iii) I a„ ) is strictly positive and strictly increasing, if a, > 0. 



Niels B. Kammerer and Wolfgang Stummer/ Distances of branching processes 7 

Due to the linear interrelation (19), the monotonicity carries over to the sequence (bn ) (p € [0, oof, 

V / n<=N 

q G]0,oo[) in the following way: 

(p-iv) b { y ] = -a x < 0, 
(p-v) bT 9 ' = p - a\, if q = p\, 



(p-vi) (bn' ) is strictly decreasing, if q < /3a, 

(p-vii) (bn 9 ) is strictly increasing, if q>j3\. 



Notice that the sign of &„ might not be same as the sign of a 9 (see e.g. (p-i), (p-iv)). Finally, for the 
remaining case one trivially gets 

(p-viii) al 0) = -/3 A , 6r ( f' 0) = e-f>* -p-ctx (p > 0). 

Moreover, for (/3„4, /3-n, a, a, ot-Hi A) € P X ]Q, 1[ an d q G]0, oo[ we shall sometimes use the function 

S, {9 \x) := q-e'-fo, x e M, (23) 

which has the following obvious properties: 

(p-ix) £^ 9j is strictly increasing, strictly conxex and smooth, 
(p-x) linxr-^-oo £ x 9 \x) = -/3a < 0, lim^oo C x 9 \x) = oo. 

With these auxilliary basic facts in hand, let us now start our detailed investigations of the time-behaviour 
n n- Hx(PA,n\\PH,n) f° r the exactly treatable case (a) in Theorem 2.2. 

3.1. Detailed analysis of the exact values 

(aNI) The non- immigration case ((3a, PH, a A, a n,^) G T'nixJO, 1[: 

Recall that for this set-up we derived q x := Pa^h ^ ^ ana - Pa : ~ a A a ]< = 0. According to Lemma 
A.l, one has q x < fix and thus, (a 9x ) is strictly negative as well as strictly decreasing. Furthermore, 



because of (p-ix), (p-x) and a 9x < 0, the function £ 9x hits on ] — oo, 0] the straight line id(x) := x once and 
only once. Consequently, ( a 9x ) converges to the unique solution Xq x e] — /3a, o-i [ of the equation 

Sx^Hx) = Qx -e x -px=x, x<0. (24) 

Summing up, we have shown the following detailed behaviour of Hellinger integrals: 
Proposition 3.1. For all (Pa, /3-h, a_4, a.%, A) £ Pjv/ x ]0, 1[ and all initial population sizes luq € N £/iere /io/ds 

(a) flA(iU,i||iW) = cxp{(/3^/3^ A -A/3^-(l-A)/3„)x } < 1, 

(b) the sequence (Hx(PA,n\\P-H,n)) nm given by 

Hx(PA,n\\P-H,n) = exp|a„ A) wo| =: Vx, n 
is strictly decreasing, 

(c) lim Hx(P A ,n\\PH,n) = exp{4 9f) w }e]0,l[, 

(d) lim - log#A(P4,nlI-PH,n) = . 

ra— >oo n 
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(aEF) The "equal-fraction-case" (fj A , [3 n ,a A , a n , A) € "Psp.i x]0, 1[: 

Again, one has qf := /3 A f3aT~ > 0. Furthermore, pf := a\a\[ > 0, which leads to the abovementioned 

/ E E\ / E\ / I E\\ 

relation 6„ ' = %^ On x . Hence, the results about the sequence [an ) coincide with those of the 

non-immigration case. This implies also that the sequence (^fe=i a fc ) * s strictly negative, strictly 

decreasing and converges to — oo. Hence, we get 

Proposition 3.2. For all (/3 A , f3-H,a A ,a-u, A) £ "Pgp.ixjO, 1[ and all initial population sizes luq € N £/iere 
/jo/ds 

(a) H X (P A1 \\P HA ) = exp{(/3 A /3^ A - A/3.4 -(1-A)/3h) L, + |£) } < 1, 

(6) the sequence (Hx(P A , n \\pH,n)) nen gwen 6y 

ffA^Anll-PK.n) = cxp J o^ ) u Q + |± £ a^ > 1 =: 7 A ,„ 

is strictly decreasing, 

(c) lim FA(-P4,nll P -H, n ) = , 

71— >00 

(d) lim - \ogH x (P A . n \\P n . n ) = %* xf ] . 

n-s-oo n p A 

Remark 3.3. For the (to our context) incompatible setup of GWI with Poisson offspring but nonstochastic 
immigration of constant value 1, an "analogue" of part (d) of Proposition 3.2 was established in Linkov and 
Lunyova [52]. 

3.2. Detailed analysis of the lower bounds 

In this section we assume (f3 A , f3 H ,a A ,a H7 \) <G CPsp\"Psp.i)x]0, 1[ and thus a A > 0, a H > 0, f£ ^ |^, 
7^0, jU(a;) > 0, fu{x) >0(i£ [0, oof). Concerning (15), let us derive a lower linear bound </>^(-) of <j>\{-) 
which is optimal. In order to achieve this, one can use the following straightforward properties of <p\{x), 
x€ [0,oo[ (cf. (10)): 

(p-xi) <fr\(0) = a\a^[ — a\ < (cf. Lemma A.l), with equality iff a A = an (together with j3 A =^ fiu)- 
(p-xii) <f/ x (x) = Xf3 A (f A (x)f- 1 (f n (x)) 1 - X + (l-X)l3n(fA(x)) X (M(x)r X -px > -ft. 
(p-xiii) hm x _^oo <t>\( x ) = P A Ph ~Px < (cf. Lemma A.l), with equality iff /3 A = /3-h (together with a A ^ ecu). 
(p-xiv) <t>\(x) = — A(l — A) (f A (x)) ~ (/h(x)) _ ~ 7 2 < 0, i.e. the function 4>\() is strictly concave; no- 
tice that 4>'x(0) = Xf3 A (a A /ccn) ~ + (1 — A) /3-h (a A /a-n) — fix can be either negative (e.g. for 
(Pa,0H,<xa,<*H,X) = (4,2,3,1,0.5)), or zero (e.g. for (f3 A , f3 n ,a A , a H , A) = (4,2,4,1,0.5)), or positive 
(e.g. for (j3 A , /3-h, a A , a-u, A) = (4,2,5,1,0.5)). Accordingly, the strict decreasingness and continuity of 
4>'\{') as well as (p-xiii) imply that <j>\{-) can be either strictly decreasing, or can obtain its global 
maximum on ]0, oo[, or - only in the case f3 A = /3-h - can be strictly increasing. 

(p-xv) lim^oo \6\(x)-(f\ + s\x)) = for fx := \a A \j£) -1 + (1-A)a w (|£J 

ft A fiu ~ flx < 0; notice that sa = iff/3.4 = /3h (together with a.4 ^ an)- Furthermore, 0a (0) < f\ (cf. 
Lemma A.l). If a A = a-u (and thus fj A ^ /3-h) then the intercept f\ is strictly positive, whereas for the 
case a A ^ a-u the intercept 7"a can take any sign (take e.g. (f3 A , /3-h, a A , a-u, A) = (3.7, 0.9, 2.0, 1.0, 0.5) 
for fx > 0, (J3 A , fa, a A , a H , A) = (3.6, 0.9, 2.0, 1.0, 0.5) for f\ = 0, (p A , *3 n , a A , a H , A) = (3.5, 0.9, 2.0, 1.0, 0.5) 
for r x < 0). 

From (p-xi) to (p-xv) it is easy to see that for all current parameter constellations the particular choices 

p{ := a A a]^ x > 0, q{ := p A p n ~ X > 

- which correspond to the choices 

r L x := a A a n - x - a x < 0, s{ := [3 A /3^ A - /3 A < 



^ -1 



and sa := 
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in (15) (and at least one of the two last inequalities is strict) - lead to the tightest lower bound B^ n for 
H\{PA,n\\pH,n) hi (21)- This situation coincides partially with those in Section 3.1. Formally, p\ = pf and 
<1\ = 1\ i Du t because of 7 7^ the relation b„ A ' 9a = %£ a„ A is in general not valid anymore and has to be 
replaced by the relation (cf. (19)) 

Hence, for a better distinguishability and easier reference we stick to the L— notation here. Nevertheless, 
the behaviour of the sequence I a n q A ) coincides exactly with that of the sequence ( a n 9x ) in the 

V / tiGN V / TieN 

Subsections 3.1(aNI), (aEF). In particular a„ ) is strictly negative, strictly decreasing and converges 

(a L ) (a L ) 

to the unique solution Xq x £] — 00, a\ [ of the equation 

^(x) = <&-f?-Px = x, x<0. (26) 

Consequently, because of (25) and }f^ x,qx ' — a\a\^ — a\ < (cf. (18)), the sequence 



(tf*^) 



is strictly negative and strictly decreasing. As in Subsection 3.1(aEF), we obtain 

raGN\{l} 



Proposition 3.4. For all [Pa, fin, a Ai a u, A) G ('Psp\'P,sp,i)x]0, 1[ and all initial population sizes wo G N 
there holds 



(a) Bi, := cxp{(^/3^ A -A^-(l-A)/3„)^ + (a>^ A -A^-(l-A) a „)} 

(b) the sequence (Bx, n ) n&l of lower bounds for (H\(PA,n\\Pu,n)) n gs 9 iven 

n 

by B x n := cxp < a„ A ujo + > 6^ A ' A > is strictly decreasing 

(c) lim Bf n = , 



fe=i 



(d) lim - logB^„ = V -\ (xf x) +j3 x ) -a x . 

n-s-oo n '9a ' 

5.3. Detailed analysis of the upper bounds 

As above, we again assume (/3 a, Pn, otA, QM] A) G ("Psp\7'sp,i) x ]0, 1[ throughout this section. In contrast 
to the treatment of the lower bounds in Section 3.2, the finetuning of the upper bounds is more involved. 
Because of the strict concavity of the function <f>x(-) (cf. (p-xiv)), there is in general no overall best linear 
upper bound of 4>x(') within the framework (15). Different reasonable goals might lead to different reasonable 
choices of p%, q^ (and thus of r%, s^) which might imply different behaviour of the corresponding sequence 

Bj of upper bounds in (21). This can be conjectured from the following immediate monotonicity 

properties: 

(p-xvi) < q x < q 2 =3- a ( n qi) < a { n q2) for all n G N. 

(p-xvii) Trivially, b„ = b n = — ax- In contrast, let p G]0,oo[ be fixed; then, < q\ < q 2 =>■ &!f' < 

&i p ' 92) for all n G N. 
(p-xviii) Let q G [0, 00 [ be fixed. Then, < pi < p 2 =^> b^' q) < bi? 2 ' q) for all n G N. 

(p-xix) < Pl < pa, < qi < q 2 => b% uqi) < b ( r m) for all n G N. 

(p-xx) For the case < p\ < p 2 , < q 2 < q± there is in general no dominance assertion for bn , bn'' 92 ' 
which holds for all n G N; take e.g. (/3a, Pu,a. A , a H , A) = (1, 0.6, 3, 3, 0.5), pi = 3.4641, qi = 0.7785 (for 
which A / (-) corresponds to the secant line through the points ^a(O) and <fix(l)), as well as p 2 = 3.4857, 
q 2 = 0.7746 (for which <f>V (•) corresponds to the asymptote of <j>x), and inspect the first six values of of 
the corresponding b n — sequence. 

The properties (p-xvi) to (p-xx) have corresponding effects on the behaviour 

(Px> 9a ) ^ Bxn = mm 1 ex P 1 a ™ ^o + Z)3b=i "k X f ) If (°f- (21)) of the upper bounds. For instance, 
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for any fixed admissible intercept p 1 ^ one would always choose the smallest admissible q x in order to achieve 
the smallest possible upper bound; due to (p-xiv) this implies that on the ultimately relevant subdomain No 
the linear function <fi x (-) should hit <fi\(-) in at least one but at most two points (tangent or secant line). 
Furthermore, we require for the rest of the section that p\ > and q x > 0, because otherwise r x < </>a(0) 
and s 1 ^ < s\ (cf. (p-xv)) which contradicts to the nature of linear upper bounds of <j> x . 

The (only partially restricted) choice of parameters p^, q x for the upper bounds B x n can be made according 

to different, partially incompatible ("optimality-" respectively "goodness-") criteria, such as: 

(Ga) very good tightness for n > N for some fixed large N G N, or 

(Gb) for a fixed initial population size wjeN there holds B Xn < 1 for all n G N, or 

(Gc) there holds B Xn < 1 for all n G N and all u>o G N (strict improvement of the general upper 

bound (13)). 
For the sake of brevity, we investigate only goal (Gc) (with the exception of Subsection 3.3(a7) and Theorem 
6.3) which can be achieved if (and "nearly but not fully" iff) (16) holds; this can be seen from 

B% A = min{exp{ai ,?) w + &f M) },l} = minjexp {(<? - p x )u + (p- a x )} , l} 

and the properties (p-i) to (p-vii). Furthermore, (p-xiv) and (p-xv) imply that the slope s x := q x — (3\ in 
(15) should be greater or equal to the limit slope sx which leads to the restriction q x > P\P\[ X ■ Moreover, 
since s x < 0, the intercept r x := p 1 ^ — a x in (15) should be greater or equal to (j) X (0) and thus, p x > a\a\7 x . 
By comparing the above established lower and upper parameter-bounds, from Lemma A.l it follows that the 
case q x < j3\ automatically implies Pa ^ fa whereas the case p^ < ct\ leads to a a ^ oen- In consistence 
with (p-xiv), various different parameter constellations can lead to different Hellinger- integral- upper-bound 
details, which we investigate in the following. 

(al) The case Psp^ of all (componentwise) strictly positive (Pa, Pn,&A, a -H) with fa -£ Pa, 

a a = ol-h 

We have </>a(0) = (cf. (p-xi)), 4> x (0) = (cf. (p-xii)). Thus, the only admissible intercept choice is r x = 
= p x — a x = of a ' 9a (i-e. p\ = cx\ = a, > 0), and the minimal admissible slope which implies (15) for 
X G N is given by s^ = 0A(1 j:g A(O) = q" - p x = a[ q " ] < (i.e. q" = (a. + p A ) X (a. + Pu) 1 ' X - a. > 0). 
Analogously to Subsection 3.1(aNI), one can derive that (a„ ) is strictly negative, strictly decreasing, 

V / nGN 

and converges to the unique solution Xq a g] — oo, af x [ of the equation 

d 9?) (*) = q%-e x -Px = x, x<0. (27) 

Moreover, in the same manner as in Section 3.2, the sequence \bn X ' A ) is strictly negative and 

V /n£N\{l} 

strictly decreasing. This leads to 

Proposition 3.5. For all (PAiP'Hi a -Ai 0i 'Ht^) G Psp.2 x ]^t^[ an d a ^ initial population sizes w G N there 
holds 

(a) B% A := exp { (<# - /3 A ) u + {p u x - a A ) } < 1, 

(6) the sequence (B^ n ) n( _ N of upper bounds for (H x (PA, n \\Pn,n)) neN given 

n 

by B Xn := exp < al? A loo + /. &fc A ' 9a f is strictly decreasing, 



k=i 

(c) lim B% n ■--- = lim H x {P A , n \\P n ,n), 

n— t>oo ' n—^oo 

(d) lim - logB% n = % (xf } +p x ) -a x 



In contrast to Psp,2j the constellation Psp,3 of all (componentwise) strictly positive {Pa-iPut&A-i&'h) with 
a A 7^ ce-u, Pa t^ Ph an d §^ 7^ Ip- is divided into three main parts as follows: because of Lemma A.l one 
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gets on the domain ]0,oo[ the relation <fi\(x) = iff f A (x) = fni x ) iff x = x* := g w Zg' 4 > 0. Accordingly, 
(for reasons which will be explained below) we denote by Psp,3ab resp. Psp^c resp. Psp^d the subset of 
V S pi for which a R H ~V < resp. %*=%* GlO,oo[\N resp. a R n ~?i A G N; notice that the case °^y = 
can not appear within "Psp,3- For further investigations let us first divide the set "Psp.3ab x ]0, 1[ of quintuples 
(Pa, Ph, ola-, an, A) into two parts V S p-^ and V^° h : 

(a2) The case V^° a of all (p A , Pu,a A , a u , A) G Psp,3a b x]0, 1[ for which \fi A (a A /a n ) X - X + (1 - 
X)Ph (a A /a n ) X - px < holds 

From (p-xi) and (p-xii), one gets <fix(0) < and 0' A (O) < 0. For the latter, both the strict negativity as well as 
the vanishing can appear in the current parameter setup, take e.g. (Pa, Pu, a A , a-u, A) = (1.8, 0.9, 2.8, 0.7, 0.5) 
for (/>' A (0) = and (Pa,Ph, a A ,a H , A) = (1.8,0.9,2.7,0.7,0.5) for <p' x (0) < 0. In the current setup, <p x is a 
strictly negative, strictly decreasing, and - due to (p-xiv) - strictly concave function (and thus, the assumption 
°a i Z ( Q A < is superfluous here). In contrast to Subsection (al), one has the flexibility to choose the intercept 
p x from the nonempty interval [a A a^ x ,a\] and the slope q x from the nonempty interval [P A P\t X , Px], 
subject to the constraints that (px,q x ) ^ (ax, P\) and <j)\{x) < (p x — ax) + (<?x — P\) x (cf. (15)). Of course, 
one way to obtain a reasonable choice of intercept and slope is the search for the optimum 



(pU, <#) := argmin M I exp {e#> Wo + £ &£»*>} I (28) 

subject to the abovementioned constraints. However, the corresponding result generally depends on the choice 
of the initial population size ujq and the observation horizon n. Hence, there is in general no overall optimal 
choice of p^, q x (without the incorporation of further goal-dependent constraints such as lim.n_j.oo B Xn = 
in case of lim„_>. 00 Hx(P An \\pH.n) = 0). By the way, due to the recursive nature of the sequences in (28) 
and the nontriviality of the constraints, this optimization problem seems to be not straightforward to solve, 
in general. 

Inspired from Subsection (al), a more pragmatic but yet reasonable choice is the following: take any intercept 
P\ e [a>",«A] such that (p%-ax)+2(<t>x(l)-(p%-a x )) > </>a(2) (i.e. 2 (a A + Pa) X (a n + Puf~ X -p l l + 
ax > (cm + 2Pa) X (an + 2/3«) 1 " A ) and <# := A (1) - (p u x - a x ) + Px = (a A + P A ) X (a n + Ph)'~ X - P U X , 
which corresponds to a linear function <^> A which is 

(a) nonpositive on No and strictly negative on N, 

(b) larger than or equal to 4>x on No, strictly larger than <fix on N\{1, 2}, and equal to (j>\ at 
the point x = 1 ("discrete tangent or secant line through x = 1"). 

One can easily see that (due to the restriction (14)) not all p A G [a^a^~ , ax] might qualify for the current 

purpose. For the particular choice p x = a A a]^ and q x — (a a + Pa) (an + Ph) ~ a A a H~ one obtains 

r x = Px - "a = b { j p "' q " ) < (cf. Lemma A.l) and s^ = q^ - Px = (j>x(l) - 4>\ (0) = a[ q " ] < (secant line 

through 4>x(0) and 0a(1))- Hence, analogously to Subsection (al) one can derive that (an J is strictly 

V / n£N 

negative, strictly decreasing, and converges to the unique solution Xq a g] — oo, a^ [ of equation (27). More- 
bn x ' qx ) is strictly negative and strictly decreasing. Thus, all the assertions (a), 

/nGN\{l} 

(b), (c), (d) of Proposition 3.5 hold for all (Pa, Ph, a a, ay,. A) G ^splia ana - au initial population sizes luq G N. 

(a3) The case Vsp£ b of all (P A ,p H ,a A ,aH,\) G "Psp,3abx]0, 1[ for which X(3 A (a^/a«) A_1 + (1 - 
X)Ph (a A /a H ) X -Px>0 holds 

In this situation (which appears e.g. for (Pa, Ph, a A, a-u, A) = (1.8, 0.9, 2.9, 0.7, 0.5)) one gets from (p-xi) and 
(p-xii) the two inequalities 4>x (0) < and <^(0) > 0. Furthermore, in accordance with the arguments in the 
forefront of Subsection (a2), <f>x is a strictly negative, strictly concave, hump-shaped (cf. (p-xiii)) function. 
One can proceed similarly to (a2). Indeed, let x max := argmax^io ^(px^) which is the unique solution of 



AAa 
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A 






+ (1 - X)Pu 



( !a{x) 



= , xg]0,oo[ , (29) 



(cf. (p-xii), (p-xiv)); notice that x* formally satisfies the equation (29) but does not qualify because of the 
current restriction x* < 0. 

Let us first inspect the case ^(L^maxJ ) > ^(L^maxJ +l)j where \_x\ denotes the integer part of x. Consider the 
subcase </>A(L a -maxJ) + L^maxJ {<j>\ ( L^maxJ ) — <fr\( L^maxJ + 1)) < 0, which means that the secant line through 
^A ( L^maxJ ) and 0A(|_ x maxJ + 1) possesses a non-positive intercept. In this situation it is reasonable to choose 

i U U \ 

as intercept any P^-a X = b^ x ' qx> = r% G [4>\ ( [x max \ ) , (f)\ ( |x m axj ) + L^maxJ OMbmaxJ) - "MbmaxJ + 1))]: 

and as corresponding s/ope g^ — q:a = Oj A = s a ~ 

(\x™"*\)-q x < 0. A larger intercept would lead to a linear function </>^ for which (15) is not valid at 

L^maxJ T 1 • 

In the other subcase 4>\{[x mlix \) + a; max (<f>\ ( [x max \ ) - <j>\ ( Lx max J +1)) > 0, one can choose any inter- 
cept p^ — a\ = fef 3 A ,?a = r^ G [0A(L a; maxJ),O] and as corresponding slope q% — a\ = a^ x = s^ = 
fli lVo' il — ^ (notice that the corresponding line </>^ is on ] [iEmaxJ , °o [ strictly larger than the secant 
line through 0a ( L^maxJ ) and ^(L^maxJ + 1)). 

If </>a ( L^max j ) < 0A(L a; maxJ +1), one can proceed as above by substituting the crucial pair of points 
(L x maxJ, Ij^maxJ + 1) with ( [^maxj + 1, L^maxJ + 2) and examining the analogous two subcases. 

With the accordingly derived pV, q^ one gets in all four (sub)cases exactly the same kind of behaviour of the 
sequences ( a„ ) , (&n ' ) as in Subsection (a2). Hence, all the assertions (a), (b), (c), (d) of 

V / nEN V /n6N\{l} 

Proposition 3.5 hold for all (/3a, /3%, aj\,, O-Hi A) G V S p 3b and all initial population sizes wo G N. 
(a4) The case T^sp^c of all (componentwise) strictly positive (f3 A , Pu , &A > a u ) with a A ^ a-n, 

Pa ^Pu,^^^ and f^ti e]o,oo[\N 

The only difference to Subsection (a3) is that the maximum value of </>a(0 now achieves 0, at the positive 
non-integer point x max = x* — g W Zg A g]0,oo[\N (take e.g. (/3j., /3-n, a A , a-n, A) = (1.8,0.9, 1.1,3.0,0.5) as an 

example). Due to (p-xi), (p-xii) and (p-xiv) one gets automatically A/3.4 (au/a«) ~~ + (1 — A)/3^ (a^/a-n) — 
j3\ > for all A G]0, 1[. This situation can be treated exactly as in (a3). Consequently, all the assertions (a), 
(b), (c), (d) of Proposition 3.5 hold for all (/3^, fin, ot-Ai a u, A) G / Psp,3c x ]0, 1[ and all initial population sizes 

uq G N. 

(a5) The case "Psp,3d of all (componentwise) strictly positive {PaiPui&Ai&h) with a A ^ a-u, 
Pa^ 0u, — ^ ir- and T~ a R A e N 

The only difference to Subsection (a4) is that the maximum value ot (f>\(-) now achieves at the integer point 
a; ma x = x* = f^ff 4 G N (take e.g. (fa, fin, a^., an, A) = (1.8,0.9,1.2,3.0,0.5) as an example). Under the 

restriction that exp < a„ A ujq + Y^k=i &fc ' > < 1 for all n G N and all wo G N, our method leads to the 

choices r^ = as well as s^ = 0. Consequently, BV = 1, which coincides with the general upper bound 

(13), but violates the abovementioned desired goal (Gc). 

However, by using a conceptually different method we can nevertheless prove the convergence 

lim H x (P A . n \\P H . n ) = (30) 

n— >oo 

(which will be used for the study of entire separation below). This will be done in Appendix A.l. 
As a next step, let us investigate the last possible parameter constellation: 
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(a6) The case Psp^ of all (componentwise) strictly positive (/3a, fitti &A> ^"H) with a. a ^ a-Ht 

This is the only case where ^a(') is strictly negative and strictly increasing, with lim^^oo 4>\(x) = lim 2: _i. 00 (j)'\( x ) 

0, leading to the choices r^ = as well as s^ = under the restriction that exp < a„ u>o + Sfe=i ^fc ?A f — 

1 for all n £ N and all ojq £ N. Consequently, B^ n = 1, which is consistent with the general upper bound 
(13), but violates the abovementioned desired Goal (Gc). Unfortunately, the proof method of (30) can't be 
carried over to the current setup (see Appendix A.l). 

(a7) Alternative bounds for :Psp,2 U Psp,3ab U Psp^c U Psp.ad 

Within this last subsection, let us exceptionally ignore the Goal (Gc). Correspondingly, for the derivation 
of an upper bound B^ n one can use the asymptote of ip\ given in (p-xv) to end up with p^ := f\ + a\ = 

Xa A (f^)^ 1 + (1 - A)a« (|^) A as well as q% = s~ x + /3 A = P A ^ X . Clearly, p" > p L x = a A a^ X by 
Lemma A.l and q^ = q^. Furthermore, q^ < j3x and thus (16) holds, since we have excluded "Psp.4- However 
- depending on the choice of (PAi^Hi a Ai a 'H) _ the intercept f\ = p^ — a\ may become strictly positive, 

and hence 

Bl x := exp {a^ } co + bf> * >} = exp { (^ - &) • ca + p" - a x } 

may become larger than 1. However, according to properties (p-ii) and (p-vi) the sequence 

„-► 2?£ n := expjaF^o + ^bf'^^cxpjal^ + l^^ap+^^-^j.n} 
may become smaller than 1. Let us therefore define for all n £ N and all A <s]0, 1[ 



{ B ln^} 



which can be used as an upper bound for the case "Psp,3<ix]0, 1[. 

For the other cases ("Psp,2x]0, 1 D U "P S p- 3a U 7 5 S p > 3b U ('Psp,3cX]0, 1[) all the assertions (a),(b),(c) of Proposition 
3.5 remain valid for replacing B^ n by the improved upper bound 

<f Pr ■■= min {*£„,£)£} < 1. 

In fact, for all these parameter classes there are concrete examples such that the upper bound S A '™ pl really 

improves the upper bound B^ n for alln £ N (i.e. B^ n < B% n ). For'Psp,2x]0, 1[ take e.g. (J3a, 0H, &A, &H> A) = 
(0.8,0.6,2,2,0.5) and u = 10, with p% = 2.021, iff = 0.693, instead of the proposed choice p v x = 2 
and q l l = 0.698. For ^g P f 3 ° take e.g. {0a,Ph,(XA,<xh,X) = (0.8,0.6,2,1.9,0.5) and uj q = 10, with pf = 
1.963, qU = 0.693, instead of the proposed choice p^ 1 = 1.949 and q^ = 0.696. For 'P S p > 3b take e.g. 
{^A>^H, a A, a n, A) = (0.8,0.6,2, 1.1,0.5) and wo = 10, withp^ 7 = 1.501, q^ = 0.693, instead of the (amongst 
others proposed) choice p^ = 1.483 and q x = 0.699. For Psp^cXJO, 1[ take e.g. {PaiPu, &A> a ~Hi A) = 
(1,1.5,2,1.8,0.5) and ujq = 10, with p^ = 1.960, q x = 1.225, instead of the (amongst others proposed) 
choice p% = 1.897 and q% = 1.249. 

3-4- Asymptotic distinguishability 

For each n £ No, let (Cl n , F n ) be a measurable space equipped with two probability measures P n , P n . The 
following two general types of asymptotic distinguishability are well known (see e.g. LeCam [42], Liese and 
Vajda [47], Jacod and Shiryaev [29], Linkov [51], and the references therein): 

(CEa) the sequence (P„)„gN is contiguous to the sequence (P n ) ng N - hi symbols, (P n ) < (P n )) - if for all 
sequences A n £ T n with limn^oo P n (A n ) = there holds lim„_ 5 . 00 P n (A n ) = 0. 
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(CEb) the sequences (P n )neN an d (-Fn)raeNo are entirely separated (completely asymptotically separable) - in 
symbols, (P n ) A (P n ) ~ if there exist a sequence n m f oo as m f oo and for each m £ No an A Um £ Tn m 
such that limm^oo P„ m (A„ m ) = 1 and lim m ^oo P Um (A nm ) = 0. 

The corresponding negations will be denoted by <3 and A. As demonstrated in the abovementioned references 
for a general context, 

(CEb) holds iff liminf„_ ! . 00 H\ \.P n \\Pn ) = for some (or equivalently, all) A <G]0, 1[; furthermore, 
(CEa) holds iff liminf An {liminf^oo H x (K\\K)} = 1. 

Combining these results with the respective part (c) of Propositions 3.1, 3.2 and 3.5 as well as the connected 
investigations of Subsections 3.3(a2) to (a5), we obtain the following 

Corollary 3.6. (a) For all {PaiPhi&AiOi-h) € (V S p,i U V S p.2 U V S p,3ab U V S p,s c U Vsp,3d) and all initial 
population sizes loq £ N, the corresponding sequences {PA,n)neN an d (Py..n)nefio are entirely separated. 

(b) For all (PaiP'Hi&Ai&'h) €E Vni with /3^ < 1 and all initial population sizes luq £ N, the sequence 
(PA,7i)neN is contiguous to {Pu,n)neN - 

(c) For all ((3a, Phi a A, &■%) € Vni with 0a > 1 an d all initial population sizes luq £ N, the sequence 
(-P.A,n)neNo * s ne ither contiguous to nor entirely separated to (P-H in ) ng N . 

Remarks 3.7. (i) Assertion (c) of Corollary 3.6 contrasts the case of Gaussian processes with independent 
increments where one gets either entire separation or mutual contiguity (see e.g. Liese and Vajda [47]). 
(ii) By putting Corollary 3.6(b) and (c) together, we obtain for different "criticality pairs" in the non- 
immigration case Vni the following asymptotic distinguishability types: 

(Pa.u) <>(p«,„) ifpA<\,Pu<i; (PA,n)<>(Pu,n) if0A<i,0H>i; 

(PA,n)<t> (P-H,n) if Pa > E Pn < 1/ (*V) <> Pn, n) and (P A ,n) A(P W ,„) if P A > 1, Pn > ^ 

in particular, for Vni the sequences (-P.A,n)neNo and (Pn,n)neN are not completely asymptotically inseparable 

(indistinguishable) . 

(Hi) In the light of the abovementioned (CEa) resp. (CEb) characteriztions by means of Hellinger integral 

limits, the finite-time-horizon results on Hellinger integrals given in Theorem 2.2, Section 3 and also in the 

following Section 4 can loosely be interpreted as "finite-sample (rather than asympotic) distinguishability" 

assertions. 

4. Closed-form bounds 

Depending on the parameter constellation, we have given bounds respectively exact values for the Hellinger 
integrals, which can be obtained with the help of recursions (17) (together with (19) respectively (p-viii)) 
which are "stepwise fully evaluable" but generally seem not to admit a closed-form representation in the 
observation horizons n; consequently, the exact time-behaviour of (the bounds of) the Hellinger integrals can 
generally not be seen explicitly. To avoid this intransparency (at the expense of losing some precision) one can 
approximate (17) by a recursion that allows for a closed- form representation. Accordingly, we shall employ 
(context-adapted) linear inhomogeneous difference equations 

S := ; a n := £(a„_i) + p n _i, n £ N, with (31) 

£(x) := c + d-x, a; e]- oo,0], (32) 

p„_! := K x ■ x™- 1 + K 2 ■ i/"" 1 , n £ N, (33) 

for some constants c £] — oo,0[, d £]0,1[, K\ , K 2 , x, v £ R with < v < x < d. As usual, one gets the 
closed-form representation 

1 J1TL slTh Tt sJTl Tt 

a n = a h n om + c n with a*r n = c ■ ^—r and c n = K x ■ - ~ - + K 2 ■ ~ V (34) 

1 — a a — x a — v 

which immediately leads for all n £ N to 

^ ak \d-x ■ d-v l-d/' l-d (d-x)(l-x) (d_ „)(!_„) l-d' [ ' 
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Notice that for the special case K 2 = —K\ > one has from (33) for all integers n > 2 the relation p„_i < 
and thus a n — aJj om < 0, leading to 

c„ < and YJ c„ < . (36) 

fe=i 
In the following, we appropriately apply (31)-(35) to the different parameter contexts of Section 3. 

4-1- Closed- form lower bounds 

Let (fa-, fa,a^, ct-u, A) € Vx]0, 1[. We have seen in the Sections 3.1 and 3.2 that the determination of the 
exact values and the lower bounds had (more or less) identical structure: choose q* := q^ = qf = f3\fi]^ X > 

0, compute the sequence ( a„ A I by the nonlinear recursion (cf. (17), (23)) 

V / TieNo 

aS^r-0 ; af x) := ^(ajg) , ^N, (37) 

choose p* '■= P\ = pf = c^A a u — ^' compute (cf. (19)) 

and finally end up with (cf. (21), (20)) exp I an x ujq + J2k=i "k A ?A f wmcn i s either interpreted as bound 

B^ n in the parameter case (fa, fa, ax, a-H, A) G ('Psp\^ , sp.i) x ]0, 1[ or as exact value V\. n in the parameter 
case (fa,fa,a^,a-n,X) £ ('PniU'Psp.i) x ]0, 1[ (where we achieved some further simplifications above). Since 



is strictly negative, strictly decreasing and converges to the unique solution x A e] — oo, a^ A [ = 

- oo, q* - fa[ of the equation (cf. (26), (24)) 

$\x) = qf ■ e x - fa = x, x<0, (38) 

we use the following approximative linear recursion in order to obtain a closed-form lower bound for both 
(here identically treatable) cases L, E: 

^ } :=0 ; «<*> := ^(Sl) + *&>, n€ N, (39) 



i.e. we replace the nonlinear function £ a 9a (x) = q* ■ e x — j3\ by the tangent line of £ A ?A at x = Xq x defined 

^• T (x) := xfi ] (l-qt- e 3 ^) + q* ■ e*^' • * , x e [a^.O] , (40) 

and reduce the error we face by adding the "correction-term" 

d§ -r(»? ?) ) a, (^" e4,?) ) 9B " 1>a (41) 

In other words, by means of the two functions on the domain [0, oo[ 

q » d^ T := q ■ e^ q * T® := § • e^ ■ (x^)" = ^ • (x^f (42) 

we use (31), (32), (33) with constants d := d<«*)- T e }0, 1[, c := 4 9a) • (l - d^^'A e]-oo,0[, K x :=T^t) > 

0, >t := [dy l \>> T \ <E R\{d, 1}, i^2 := 0, ^ := 0. Let us first present some fundamental properties which 
will be proved in Appendix A. 2: 

Lemma 4.1. For all (fa, fa, aj(,au, A) <G ?x]0,l[ there holds: 

(a) al?* ) < a [ f x) , forallneN. 

(b) The sequence I a n A j is strictly decreasing. 

V / nGN 

(c) lim al" A) = lim ai 9 * 5 = x^^- 
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Applying Theorem 2.2. Lemma 4.1 as well as the formulae (19), (34) and (35), one gets 

Theorem 4.2. For all (Pa,P'H, a A, a«, A) G ?x]0,l[ and all initial population sizes loq G N £/ie following 
assertions hold: 

(a) for all observation horizons n G N the Hellinger integral can be bounded from below by the closed-form 
bounds Hx(PA,n\\Pn,n) > C^ n given by 



n L 



exp U 9t) 



w 



Pt 



d(s*)' T 



1 



d dt), 



g* 1 - d(9*)- T 
+ C^* } • w + J?i 9 * } I , where for all n G N 






«A 



C 



(9*) 



r(«*) 



»(»?) 



1 - d(9x) 



'it 



- d^ 




l + d(«)> 



> o. 



(43) 
(44) 
(45) 



(fy) i/ie sequence I Cjr ) is strictly decreasing. 

(c,) /or aH observation horizons n£N 



C* 



A,n 



< 



B£ B) i/ (/3^,/3«,a^,a-H,A) G (Psp\^sp,j)x]0, 1[, 
Vx,», tf {Pa, faction, X) G (P ffi U%i)x]0,l[. 



W 



0, b */(/3^,/3„,a^,a«,A)e(7' W /) c x]0,l[, 

exp |wo x 9a } } > °> */ (Pa,Ph,(* a ,ciih,X) €VniX]0, 1[, 
w/uc/i coincides with linv^^ i?A(-P4,n||-pH,n) /or a// (/^U, /?«, a a, o-h, A) G CP\Pspu) x ]°; ![• 



lim C A .„ = 

n— >oo 



fej lim ±logC£ n = % 



M 



x a + /3a I — C*A i which coincides with 



linij^oc i log H A (fU,n | l-P-H.n) /or all (Pa, Pu, a a, an, A) £ (PniWsp.i) x]0, 1[ respectively wiift lim^oo i logBf„ 
/or aH (#A,frt,au,aK,A) € (?sp,2 U P s? , 3a6 U ? s? „, c ) x]0,l[. 

Remark 4.3. Notice that the formula (43) simplifies in the parameter case ([3 a, Pu, ola, ecu, A) € 'Psp.ixjO, 1[, 
since t/ien ii /iota's P* /q* = §^ = fp 1 anc ^ therewith (p*/q*) -/?a — q a = 0; /or i/ie case (/3.a, (3-u, a^, a%, A) G 
'PiV/xjO, 1[, one can euen use i/ie stronger relation p x = = a\. 

(a*) 

In order to get an "explicit" lower bound which does not rely on the implicitly given fixed point x A , one can 

replace the latter by a close explicit lower approximate x Q A < x Q A and proceed completely analogously, 
leading to a smaller lower bound (say) C_ x n < C Xn in assertion (a) of Theorem 4.2; in the corresponding 



assertions (b), (c) and (d) one then has to replace C x n by C_ x n and Xq a by Xq a 



choose 



where 



Sit) 






l -*t )-\ (1-«*J -2-9* • 



=h(9? 



(«*-/9a) 



For instance, one could 
(46) 



Kqt) 



f max{-/? A ; ^^} , if 9* < 1, 

I -/?a, $qt > 1; 

this will be used as an auxiliary tool for the diffusion-limit-concerning proof of Lemma A. 3(c) in the appendix. 

If q* < 1, the ter 



m % — sr represents the existing negative intersection of the tangent of £ A at x = and 



-'/a 



Sit) 



Sit) s„ 



.(«*), 



the bisectrix. Clearly, h(q*) < Xq x . By (46), x A is the unique negative solution of Q x A (a;) = x with the 
quadratic function 
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Notice that Q { ^\o) = £>\o), %- (0) = %- (0), ^- (*) < ^- (») for all x G [4* x) ,0], and 
thus Q ( f\x) < $*\x) for all x G [x£*\o[, which leads to the desired x^ 9x) < x£*\ 

4-2. Closed-form upper bounds 

In order to achieve closed- form upper bounds, we principially proceed as in the previous Section 4.1. How- 
ever, the situation is now more diverse since we have to start from Section 3.3 which carries much more 
"nonuniqueness" respectively variety than the corresponding Sections 3.1 and 3.2 which we used as a starting 
point for the investigations in Section 4.1. 

Notice first that for the subcases "Psp,3abx]0, 1[ and "Psp^cXjO, 1[ (cf. Subsections 3.3(a2),(a3),(a4)) one can 
achieve a closed-form upper bound without further investigations: if one chooses q x = f3\ (and thus, the slope 
s^ =0), then by properties (p-i), (p-v) one has a n = (i.e. recursion (17) is trivial), &„ = p — a\ < 
and hence n 

B£„ = expja^o + ^6^} = exp {n ■ (p u x - a x )} ^ . (47) 

fc=i 
However, there might exist (and for 'Pgp~3 a definitely exists) choices (p x , q x ) which lead to (fully or eventually 
partially) tighter upper bounds B^ n but for which the non-linear recursion (17) is nontrivial. Such potential 
cases, for which in particular < q x < j3 x and < p 1 ^ < a x holds, will be treated in the following; since 
the parameter constellation (/3^, /3-%, a A , a-n, A) G ("Psp,3d U T^sp.-OxJOi 1[ does not meet this requirement, let 
us fix (f3 A ,f3 n ,a A ,ay,, A) G (7>ni U 7>sp,i U Psp.2 U Psp,3ab U Psp,3c)x]0, 1[ where we also include the two 
setups Pni U "Psp.i f° r which we want to replace the recursive, non-closed-form exact values by closed-form 
upper bounds. For this situation, we determined recursive upper bounds respectively exact values in a (more 
or less) identical structure which is also very close to the one given by (37) to (38): choose q x for G = U 

respectively G = E subject to the corresponding parameter case which leads to a^ = s x = q x — /3\ < 0, 
compute the (rest of the) sequence ( a„ x ) by the nonlinear recursion (cf. (17), (23)) 

V / n£Nn 



neN 

^ G) :=0 ; a£* 



i q?) {ai q !l), neN, (48) 



choose p x subject to the corresponding parameter case and evaluate 

$U?) = tf a (f.) + fi Px _ a neN; (cf . (19)) 

which leads to the desired bound B Xn = exp < a„ ujq + YHk=i "k A ' ?A f ( c ^ P ar ^ ( D ) °^ Proposition 3.5). 
According to (p-ii) , the fundamentally important sequence ( o^ ) is strictly negative, strictly decreasing 

V / riGN 

and converges to the unique solution Xq x g] — oo, af x [ of the equation 

tf\x) = q^-e x -l3 x = x, x<0. (cf. (27), (24)) 

For an upper bound of the sequence a n x we introduce the recursion 
,(??) ._ n . ^ 



a n := U ; a r , 



if^SX) +P& nen, (49) 



i.e. we replace the nonlinear function £ x 9x (x) = q x ■ e x — fi x by the secant line of £ x across its arguments 

/ G\ 

Xq X and 0, defined by 



„G 



(9a) („G 



^> S (*) := tf-fc + *>* -fa-™ -*, xe[ X f\0] , (50) 
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and reduce the error we face by adding the "correction-term" 



Pn-l — 



ii^T-i 



i - 



' 4'>>-( g g-ao 



< 0. 



In other words, by means of (42) and the function on the domain [0,oo[ 



,.<"> 



q ^ 



d ( q ),S ._ x ~ (g - P\) 



(?) 
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(51) 



(52) 



we use (31), (32), (33) with the constants d := d^' S e]d^< T , 1[, c := q^-Px e]-oo, 0[, Kx := -T^ < 0, 
k := d(9?)- T e]0, d[, K 2 := -K x , v := d( ? ?)' T • d^?)' 5 e]0, x[. The following fundamental properties will be 
proved in the appendix: 

Lemma 4.4. For all (Pa,Pu, a A ,a n , A) € P\{7>SJ> f MU7>sJ>4}x]0,l[ it holds 



(«?) 



(a) a„ > a,n x ' , for all n£N, with equality iff n = 1 . 

(b) The sequence ( a,\ A I is strictly decreasing. 

V / tiGN 



(c) lim a 



^ = lim a£ ?) 



,(«?) 



Applying Theorem 2.2, Lemma 4.4 as well as the formulae (19), (34) and (35), one obtains 

Theorem 4.5. For all ((3a, Ph, cm, ««, A) € (Pni U Psp.i U Psp^ U Vsp,sab U Psp,.?c)x]0, 1[ and aZZ initial 
population sizes ujq € N i/ie following assertions hold: 

(a) for all observation horizons n £ N i/ie Hellinger integral can be bounded from above by the closed-form 
bounds Hx(PA,n\\Pn,n) < C\ n given by 



C\, n ■= exp ^ z 




W0 ^l-d(€).« 


(l-(dta?^) n ) + 


"4 (a- 

L<?A v 


f4^)- 


"A 


• n 




— C, n x ■ loq — d n A > , where for all n G N 




(53) 


7(9?) ._ r (??) 

Sri ' x 




> 


and 


(54) 


3?) r(^)-pG. d ( ? ?),T 
g?(l-<K>' T ) 


"l - ( d foS?),s d ta?).T) B ! _ (>?)■*)" 1 - ( d (??)^)" - 


>0 . 


(55) 



1? 



T/ie parameters < q 1 ^ < f3\. < p^ < a\ can be chosen subject to the restrictions explained in the 
parameter- adequate Subsections 3.3(al),(a2),(a3),(a4) (for G = U) respectively Subsections 3.1(aNI),(aEF) 
(for G = E). 

(b) the sequence ( C? J is strictly decreasing. 

V ' / men 

(c) for all observation horizons n£N 






> 



with equality iff n = 1. 
(d) 



lim CZ T 

n— »-oo 



B% n , if (Pa, Pu, ola, «w, A) € (V SP , 2 U 7><jP, S a6 U 7> S p, 3c )x]0, 1[, 
Vx,n, */ (/3.4,/3«,a4,a«,A)€ (? ffl UP 5 p,i)x]0,l[, 



0, tf (pA,Pn,a A ,a n ,\) &Vsp\{Vsp,sd^^sp,M\^l 

cxpjwoa^'} >°, */ (0A,Pu,oiA.,aH,X) eV NI x]0,l[, 



lim F A (PA in ||P Kin ). 



(e) lim — logC 

n— >oo 72 



G 

A, 71 



Ha. 



(4 



(<i[ 



U\ 



which coincides with 



linin^oo i log Hx(PA,n\\Pu,n) for all (Pa, Pu, a a, uuA) G (PniWPsps) x]0,l[ repectively withlim r , 
for all (p A ,Pn,a A ,an,X) G (7>sp 2 U P S P, 3a6 U 7> S p,3c) x]0,1[. 



logBF, 
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Notice that the strict positivity in (54) and (55) can be easily seen from (36). 

Remark 4.6. The formula (53) simplifies in the parameter case (,$4, /3^,a^4,a^, A) £ "Psp.i x ]0,l[ which 
results in ^i- (/3\ + Xq J — a^ = ^ • Xg . For the case (/3^, /3%, a^i, a-^, A) <G Pjv/X]0, 1[, one can even 
use the stronger relation pf = = oc\. 

(q G ) 

In order to get an "explicit" upper bound which does not rely on the implicitly given fixed point x 

(G <G {U,E}), one can replace the latter by a close explicit upper approximate Xq > Xq and proceed 

q 

completely analogously, leading to a larger upper bound (say) C x n > C Xn in assertion (a) of Theorem 4.5; 

s~i r E \ / E \ 

in the corresponding assertions (b) , (c) and (d) one then has to replace C Xn by C x n and Xq x by x 9 x . One 
possibility along these lines is the choice 



l-o G x )-J(l- q fY-2- q f-{ q f-p, 



(56) 



■" € 

which is exactly the unique negative solution to the quadratic equation 

Q ( f\x) ..= q ±. x * + q G. x + c£-p x = x. 

<q a ) ( q G ) —(q G ) 

By inspection of the first two derivatives, one gets Q x x (x) > £,\ (x) for all x € [—00, 0[, and thus x > 



x 



(q^) . . . . (a ) r* 



A . Notice that the additional fundamental requirement x < (] x — j3\ holds for parameter constellations 

/ G \ / G \ / G \ j z^i 'A / G \ 

for which x 9 x > — 1 for any Xq A < Xq x , since then one has — -^ — (x) > for all x £ [xq A , 0] . Such a 
situation will be used as an auxiliary tool for the proof of Lemma A. 3(c) in the appendix. 

5. Hellinger integral bounds in the diffusion limit 

One can show that a properly rescaled Galton- Watson process with immigration (GWI) converges weakly to 
a diffusion process X := < X t , t £ [0, 00 [ > which is the unique, strong, nonnegative - and in case of -% > ^ 
strictly positive - solution of the stochastic differential equation (SDE) of the form 



dX t = [f] - kXA alt + a\l X t dWt, ie[0,oo[, X e]0,oo[ given, (57) 



where n £ [0, oof, K £ [0,oo[, o e]0,oo[ are constants and Wt denotes a standard Brownian motion with 
respect to the underlying probability measure P; see e.g. Feller [16], Jirina [31], Lamperti [36], [37], Lindvall 
[49], [50], Grimvall [18], Jagers [30], Borovkov [7], Ethier and Kurtz [13], Durrett [12] for the non-immigration 
case corresponding to 7/ = 0, k > 0, Kawazu and Watanabe [34], Wei and Winnicki [70], Winnicki [72] for 
the immigration case corresponding to r\ ^ 0, k = 0, as well as Sriram [61] for the general case r\ £ [0, oo[, 
k £M.. Feller-type branching processes of the form (57), which are special cases of continuous state branching 
processes with immigration (see e.g. Kawazu und Watanabe [34], Li [43], as well as Dawson and Li [10] for 
imbeddings to affine processes) play for instance an important role in the modelling of the term structure 
of interest rates, cf. the seminal Cox-Ingersoll-Ross CIR model [9] and the vast follow-up literature thereof. 
Furthermore, (57) is also prominently used as (a special case of) Cox and Ross's [8] constant elasticity of 
variance CEV asset price process, as (part of) Heston's [22] stochastic asset-volatility framework, as a model 
of neuron activity (see e.g. Lansky and Lanska [39], Giorno et al. [17], Lanska et al. [41], Lansky et al. [40], 
Ditlevsen and Lansky [11], Hopfner [27], Lansky and Ditlevsen [38]), as a time-dynamic description of the 
nitrous oxide emission rate from the soil surface (see e.g. Pedersen [58]), as well as a model for the individual 
hazard rate in a survival analysis context (see e.g. Aalen and Gjessing [1]). 

Along these lines of branching-type diffusion limits, it makes sense to consider the solutions of two SDEs 
(57) with different fixed parameter sets (77, k^,<t) and (77, k-h,o~), determine for each of them a correspond- 
ing approximating GWI, investigate the Hellinger integral between the laws of these two GWI, and finally 
calculate the Hellinger integral (bounds) limit as the GWI approach their SDE solutions. Notice that for 
technicality reasons (which will be explained below) , the constants 77 and a ought to be independent of A., T-L 
in our current context. 
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In order to make the abovementioned limit procedure rigorous, it is reasonable to work with appropriate 
approximations such that in each convergence step m one faces the setup ("Pni U 'Psp,i) x ]0, 1[ (i.e. the non- 
immigration or the equal- fraction case) , where the corresponding Hellinger integral can be calculated exactly 
in a recursive way (cf. Section 3.1). Let us explain the details in the following. 

Consider a sequence of GWI (X^ m >^ with probability laws P± m > on a measurable space (fl,A), where as 

above the subscript • stands for either the hypothesis % or the alternative A. Analogously to (5), we use for 
each fixed step m G N the representation X^ := < X n m , w(fN with 

X™ := E Y n% + Y !l m \ neN, X { m) G N given, (58) 

fe=i 

where under the law P± m > 

• the collection y( m ) := i Y^{ k , n G N, k G N > consists of i.i.d. random variables which are Poisson 

distributed with parameter /3, ' > 0, 

• the collection Y( m > := < Y„ m , neN consists of i.i.d. random variables which are Poisson distributed 

with parameter a. > 0, 

• y( m ) and y(™' are independent. 

From arbitrary drift-parameters r\ G [0,oo[, k % G [0,oo[, and diffusion-term-parameter a > 0, we construct 
the offspring-distribution-parameter and the immigration-distribution parameter of the sequence ( X„ 



' nGN 

^ (3 { . m) := 1 - -£- and ai m) := pi m) ■ \. (59) 

Here and henceforth, we always assume that the approximation step m is large enough to ensure that 
/?. €]0, 1] and at least one of (yT 1 , p^' is strictly less than 1; this will be abbreviated by m G N. Let us 
point out that - as mentioned above - our choice entails the best-to-handle setup (Pni U7 : 'sp.i) x ]0, 1[ (which 
does not happen if instead of 77 one uses 7/. with 77^4 ^ Vn)- Based on the GWI X^ m \ let us construct the 
continuous-time branching process X^ n ^ := < X\ m , t G [0, oof > by 



in 



living on the state space E^ m > := — N . From (60) one can see immediately the necessity of having a to be 
independent of A, W because for the required absolute continuity in (6) both models at stake have to "live" 
on the same time-scale r t := [cr 2 m/jj . For this setup, one obtains the following convergence result: 

Theorem 5.1. Let 77 G [0, 00 [, K» G [0, oof, a g]0,oo[ and X^ m ' be as defined in (58) to (60). Furthermore, 
let us suppose that lim m _>. 00 — Xq = Xq > and denote by D([0,oo[, [0,oo[) the space of right- continuous 
functions f : [0,oo[i— > [0, oof with left limits. Then the sequence of processes [X^ m >) _ convergences in 

distribution in -D([0, oof, [0, oof) to a diffusion process X which is the unique strong, nonnegative - and in 
case of -?r > \ strictly positive - solution of the SDE 



dX t = (77 - K,X t )dt + o-yX t dW', £g[0,oo[, Xq G]0,oo[ given, (61) 

where W° denotes a standard Brownian motion with respect to the limit probability measure P». 

Notice that the condition -% > i can be interpreted in our approximation setup (59) as a, > /3, /2, 
which quantifies the intuitively reasonable indication that if the probability P,[Yn — 0] = e~ a% of having 
no immigration is small enough relative to the probability P^Y^-ik = 0] — e~^' of having no offspring 
(ttt, G N), then the limiting diffusion X never hits zero almost surely. 

The corresponding proof of Theorem 5.1 - which is outlined in Appendix A. 3 - is an adaption of the proof 
of Theorem 9.1.3 in Ethier and Kurtz [13] which deals with drift-parameters 77 = 0, k, = in the SDE 



Niels B. Kammerer and Wolfgang Stummer/ Distances of branching processes 



21 



(61) whose solution is approached on a a— independent time scale by a sequence of (critical) Galton- Watson 
processes without immigration but with general offspring distribution with mean 1 and variance a. Notice 
that due to (59) the latter is inconsistent with our Poissonian setup, but this is compensated by our chosen 
cr— dependent time scale. Furthermore, (59) is also inconsistent with the other concrete parameter choices in 
the abovementioned corresponding references. 



For each approximation step m and each observation horizon t\ . the corresponding Hellinger integrals 



Hy, P 



,(m) 



Ak 2 miJ 



P 



(m) 



U\<J 2 mt\ 



obey the results of 



(apl) the Propositions 3.1 (for r\ = 0) and 3.2 (for r\ g]0, oo[), as far as recursively computable exact values 

are concerned, 
(ap2) Theorem 4.2 and Theorem 4.5, as far as closed-form bounds are concerned; recall that the current setup 

is of type ("Pni U "Psp,i) x ]0, 1[, and thus we can use the simplifications proposed in the Remarks 4.3 

and 4.6. 



In order to obtain the desired limits of the Hellinger integrals H\ I P 



A,\_<r 2 mt\ 



\P- 



M 



U\<J 2 mt\ 



respectively of their 
closed-form bounds as m — > oo, one faces the following problems: in accordance with Section 3.1, for each 

(m) / (ttiM / ( m )\ ( m ) 



fixed 77i in (apl) one has to choose the parameters p x 

which in particular determine the fundamental sequence ( oJi ) 



'H 



1x 



en 



("1 



(cf. (17)). This enters 



M 



and the 



in the appropriate versions of part (b) in Propositions 3.1 and 3.2 respectively in form of 

't 
correspondingly arising convergences (as m . — >• oo) seem to be not (straightforwardly) tractable due to the 

recursive nature of (17). In contrast, for the closed-form bounds of Section 4 the desired convergences are 

tractable, which will be worked out in the following. To begin with, let us explicitly formulate the results 

of the application of Theorem 4.2 (where Remark 4.3 applies) and Theorem 4.5 (where Remark 4.6 applies) 

to the current setup. For this, we use the following SDE-parameter constellations (which are consistent with 

(59) in combination with our requirement to work here only on (Vni U 7-'sp,i)x]0, 1[): let Vni be the set 

of all {kAjKHiV) f° r which r\ — 0, ka G [0, oo[, k-h G [0,oo[ with ka ^ k~h', furthermore, denote by Vsp,i 

the set of all (kj(,k^,i)) for which r\ G]0,oo[, ka G [0, oo[, k-h G [0,oo[ with ka ^ k%. On Vni U Vsp.i 

there hold for m G N the useful restrictions q x G]0, 1[ and (3™ g]0,1[. For the sake of brevity, let us 



henceforth use the abbreviations a 



(m) 



„(™) 



„(m) 



(" 



(mK 



A • a^> + (1 - A) 



r (; 



r(«r ; ) = &— 



e o 



,M 






= d^ 



:= A • ^ + (1 - A) • r3%\ 



,8 _ 



>-&r-pr } ) 



and 



•' 

d (m),T ._ ( j{q x ),t _ gVn) . p x _ gy ^ e aDOve considerations, the Theorems 4.2 and 4.5 (together with 
their remarks) adapt to the current setup as follows: 

Corollary 5.2. For all (ka, K n, Vi ^) <= (Pni U / d sp,i)x]0, 1[, all t G [0, oo[, all approximation steps m G N 
and all initial population sizes Xq G N the Hellinger integral can be bounded by 



cxp < x 



(m) 




X, 



(m) V 



S 



m),T 



a 2 1 - dM< T 



d (m),T 



<7 mt 



0-2 L J 2-[a 2 mt\ ° — l<? 2 mt\ f 



(62) 



(m) 



< H,iP^ mt] 



< CXP < Xq 



X, 



p( m ) 

r U\a 2 mt\ : 

rf (m) ' S 



a 2 1 - d(™)' s 

(m) V 



1 - d (m) 



-(m) 



s \ l<? 2 mt\ 



A m) '' \rr 2 mt\ 7^' Y (m) ^ (m) L 

C ~^2 ' L CT mt \ ~ Ho- 2 mtJ A - Vy a 2 mt \ > , 



(63) 



where we define analogously to (44), (45), (54) and (55) 



»( m ) — r( m ) 



= r 
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1 - d( m ). 



>0. 



^( m ) ._ 



»M 



4. rM . ( i - ( ^:) .r.-^-- i+ r? 



1 + d(™). T 



p( m ) . 



( d M,S) n _ ( d M^) n (d (m) ' s )" 1 (l - (d(™)> T ) Tl 



>0, 



d(m),S _ a]{m),T 



1 - dW. T 



>0. 



tf (m) 



A . r( m ) 



d ll > 



i - d m )' T 



1 - ( d M^M^)™ ( d (™),5)» _ ( d (m),T)' 



1 - d{m),Sd(m),T d(m),S _ rf(m),T 



>0. 



Notice that the bounds (63) simplify significantly in the case {Kj\.,Kti,ri,X) £ 'Pjv/x]0, 1[ /or which t] = 
/io/(is. 

Let us finally present the corresponding desired limit assertions as the approximation step m tends to infinity, 
by making use of the quantities 



k\ := \ka + (1 — \)k-h > as well as Aa := a/Ak^ + (1 — A)k^ > n\ : (64) 

Theorem 5.3. Let the initial SDE-value Xq g]0, oof be arbitrary but fixed, and suppose that lim r „_>. 00 — Xq — 

Xq. Then, for all t £ [0, oof and all (kj[, K-u,r/, X) G (Vni U 'Psp i i)x]0, If the Hellinger integral limit can be 
bounded by 



Dx,t ■= exp 



Aa - k\ 



X 



V_ 

A A 



(1 - e- A ^) 4 (A A - « A ) • t 



L<£\t)-X + ^-Lf\t) 



^ ^ Hx \ P ^rnt\ 



< exp • 



(m) 

Aa - k\ 



P 



(m) 



X, 



U.\a 2 mt\ 

n 



0~ T 



(A A + k x ) 



(l_ e -|(Ax + «A)-t)_^ (AA _ KA) . t 



[/«(*)• X - 4'^fW =: ^L, 



where for all t > 



4 2) w 



(A A - k a ) _ A x -t . a _ P -Ax-t\ 



mt) 



r(2) 



^n*) 



2ct 
1 ( 


2 -Aa l 

' A A - K X 


4 1 

(A A 


y A A 
\ 2 

- KX) 


(Aa 


a 2 

\ 2 

~ KX) 



> o, 



1-e 



-A A -t\ 



> o, 



e 2 



i(A A +K A )-t _ e -AA-t e -5(A A + K A )-t /-j _ e -A A -n 



Aa 



-M, 



A A - k a 2 • A A 

^ _ g-3(3A A +K A )-t e ~ AA '* — g-1?(A A + K A H 



> 0, 



3A A + k x 



tWi 



Aa - «a 



> 0. 



(65) 



(66) 

(67) 
(68) 
(69) 

(70) 



Notice that the components L\ (t) and U\ (t) (i = \,2) do not depend on the parameter r\, and that the 
bounds (66) and (66) simplify significantly in the case (ka, K-H,r/, A) G 'Pat/xJO, If, for which r\ — holds. 



6. Power divergences and relative entropy 

All the results of the previous sections carry correspondingly over from the Hellinger integrals i/x(-||-) (A € 
]0, If) to the power divergences I\ (-||-) by virtue of the relation (cf. (1)) 
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T (P IIP 1 1-H X {P A \\P H ) 
h{PA\Pn) = A . (1 _ A) ■ 

In particular, this leads to bounds on I\ {Pa\\Ph) which are tighter than the general rudimentary bound (4) 
connected with (13). Furthermore, it is well known that in general the relative entropy defined by (2) 

I{Pa\\Ph) = Jim h(P A \\Pn) , (71) 

see e.g. Liese and Vajda [47]. Accordingly, for our context of GWI we can use (71) in combination with the 
recursive exact values respectively recursive lower bounds of Theorem 2.2 and Section 3.2 to obtain the fol- 
lowing closed-form exact values respectively closed-form upper bounds of the relative entropy I(PA.n\\Pu.n)'- 

Theorem 6.1. (a) For all (PAiP'Hi a Ai a 'H) S {PniU Psp,i), 0,11 initial population sizes luq £ N and all 
observation horizons n £ N 



I{PA,n\\Pn,n) 



3 *iMkr. 



-l)+/3w 



UJ 



^■[^■(log(|A)-l)+fa] 

3 H -\ogl3 n -l]-[^-n^ 



l-t>A 



(i - (PaT) 



^0 + ^ 



2 






(72) 



(b) For all (/3„4, f3-u, Ua, C^h) £ Vsp\Psp,i > oil initial population sizes uq £ N and all observation horizons 
neE it holds I{PA,n\\ P H.n) < E% , where 



E, 



pA-(log(^)-l)+Pn 
1-0A 



Wq 



r-PJ 



■£a 



aA-[pA-(log{^)-l)+Pn\ 
Pa(I-Pa) 



oca 



(i - {PaT) 

log(^ 



Pa 



Ha 



an 



n, ifflA^ 1 , 



(73) 



[P H -\ogi3 n -l}-[!f-n 
<*A flog i^A - Pu 



(-o + ^) 



Ot-H 



Remark 6.2. The n— behaviour of (the bounds of) the relative entropy I (Pa 
influenced by the following facts: 



if Pa = 1. 
n\\Pu,n) in Theorem 6.1 is 



(3 A ■ (log (|^) - l) + f3 n > with equality iff (3 A = Pu- 

In the case ^/lo/ (73) , there holds [ PMi-M — 

with equality iff a a = ci-h and (3 a = (3% ■ 



OCA 



LO £>\anPA) Pa 



a n > 0, 



In contrast, in order to derive (semi-) closed- form lower bounds of the relative entropy I(PA,n\\pH,n) we use 
(71) in combination with the recursive upper bounds of Theorem 2.2(b) and appropriately adapted detailed 
analyses along the lines of Section 3.3. This amounts to 

Theorem 6.3. For all (pA,PH, a A, a H) G 'PspXPsp.i, all initial population sizes ujq £ N and all observation 
horizons n £ N 



{4;f",^f c , Et h ° r }e[Q.. 



sup 

feeNo, i/e[o,oo[ 



OO 



(74) 



where for all y £ [0, oo[ we define the - possibly negatively valued - finite bound component 



jpL.tan 



(3 a In ( ^ 



+PAV \ i a 1 1 _ Qa+PaV 
+PuVj ^ iH I ccu+PuV 



1~(PaY 

i-Pa 



cj 



Pa(i-Pa) 



Pa In ( aA ti AV ) + P 



+ \a H - a a 



In 



otH+P-HV 



■Mi 



a-H+Pny 

aA+y 
otn+PnV 



H \1 
•n, 



au+Pny 



1.4 

2 



1^0 



Ml 
2 / 



+ (««-a^«)(l- = ^ 5 )-n ) 

and /or aH fc £ Nq i/ie - possibly negatively valued - finite bound component 



if Pa * 1, 
*//U = l, 



(75) 
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h k,n '■= < 



w* + 1) in (M^}j - ^( fc ) ln (£80 + ^ - f>A ■ l -^£ 
-(^(*+ 1 ) , »(feSw})-w*) ln (feS))-(*+tt) 



w 



1-&A 



f /U(fc+i) 



^( fc + 1) ln (fcM) - ^ fc ) ln (MS) + to - 



OLA) 



>U - OSH 



2 + (w + 2 f)-n] 
if Pa = I- 



(76) 



Furthermore, on Vsp,4 we set E^ ,hor := for all n £ N whereas on Vsp\(Vsp,i U Vsp,4) we defir 



E 



L.hor 



{&A + Paz*) 



lo s — r« — * I " 



+ au + Puz* 



,neN, 



(77) 



withz* :— argmax xG N I (a a + Pax) - log ( f^f^fj + 1 - («h + Pux)\ ■ In the subcases Vsp,2WPsp,3abU 



T~'sp,3c^'Psp,4 one gets even _E„ > for all wo £ N and all n £ N. In the subcase Vsp,3d, one obtains for each 

zedn 
and hence 



fixed neN and each fixed w £ N the strict positivity E% > if (-§^E^ an ) (y*) ^ 0, where y* := °£Z%* € N 



Yy E ^ an 



(78) 



(0A-M 3 . 1-(/3a)" 



W 



1-0A 



(Pa-PmV 
Pa 



aA{PA-Pu) 



(l-0A)(<XAP-H.-anPA) 



(1-/3-H) . \*A_ 



n 2 + (w + ^)-n] -(l-/^) 2 -™, 



if Pa = 1. 



Remark 6.4. Consider the exemplary parameter setup (Pa, Phi (x Ai (x 'h) = (aiai^, 1) S Psp,3d- For ini- 
tial population wo = 3 it holds y--Ej'^ " (?/*) = for all n £ N, whereas for wo 7^ 3 one obtains 

(jk E i;n m ) (I/*) ^ /or a// n £ N. 

It seems that the optimization problem in (74) admits in general only an implicitly representable solution. 
Of course, as a less tight but less involved explicit lower bound of the relative entropy I(PA,n\\Pn,n) one can 
use any term of the form max<E^* an , E^ ec , E%> hor \ (y £ [0, 00 [, k £ No), as well as the following 

Corollary 6.5. (a) For all (Pa, Pw.i a AT a ~H) £ 'Psp\'Psp,i> a ^ initial population sizes wo £ N and all 
observation horizons n £ N 

I{PA,n\\Pn, n ) > E L n > Ek := max {£&£", E^ ec , E^ hor ) £ [0,oo[ , 

with E^' hor defined by (77), with -possibly negatively valued -finite bound component E^^ a := lim„_ > . 00 E^' 1 
where 



L,ta 

y-^-oo '-'y^n > 



' *U-(i°g(f£)-i)+fa 

1-0A 



rpL,ta7i 



w 



°LA 

i-Pa 



(1 - (Pa)") 



<*A- AA-llog 



Ea 



-1 +/3w 



/3a(1-/3^) 



.n'H 



^(l-^j+owll-^ 



[/*« - log/3 w - 1] • [Sf ■ n 2 + (wo + *£) ■ n] 
+ [ .a (1 - /9«) + a« (l - ^)] ■«, 

and - possibly negatively valued - finite bound component 



n, ifPA^l, 



if Pa = 1, 



E, 



o.ii 
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i-Pa 



1-0J 



(<M + Pa) ■ log [%*ft ) - a A • log {*£ ) + Pu~Pa 
+ { T^^aI (W + M ■ log (f£fe) - a A ■ log (sa)) - ^- (i - /?„) 



-^(1 + H 



(^ + l)-log(^L)-^-log(^)+^-l].[ 
-{^ [(cm + 1) ■ log (^±L) »A ■ log (^) " 



• ^0 + 2f ■ " 2 ] 



-OA 



( 1 + ^)- 1 °g(^^))+^} 



tf#A = l. 



Por i/ie cases Vsp,i U Vsp,3ab U Tsp,3c one gets even E^ > /or aZ/ ojq G N and a^Z n G N. 
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(79) 



For the diffusion-limit of the relative entropy we obtain closed-form exact values: 

Theorem Q.Q.Within the framework of Section 5, one gets for all initial SDE-values Xq G]0, oo[, all obser- 
vation horizons t G [0, oo[ and all parameter constellations {n A , Kfi,r/) G (Vni U Psp,i) 



lim 7 P 



,(m) 



Ak 2 ™*J 



P 



(m) 



■H,l<T 2 mt] 



J^o^lM^Untj 



2 (J 2 -K T 



X -^ -(l-e-^^+r?-* , (fM>0, 



2a- 



§ • i 2 + X • t 



lim lim A P 



(„,) 



\/*l m— >oo 



A l/Ak 2 ™*] 



P 



H,\a 2 mt\ 



3 (m) 

•H,Lcr 2 mtJ 



if KA = 0, 



7. Applications 

As already mentioned in the introduction, there are numerous applications of both ingredients - power diver- 
gences resp. Hellinger integrals resp. relative entropy on the one hand and Galton- Watson branching processes 
with immigration on the other hand. In order to indicate the concrete applicability of our combinating in- 
vestigations, for the sake of brevity we confine ourselves to some issues in the context of Bayesian decision 
making BDM and Neyman-Pearson testing NPT. In BDM, we decide here between an action d-n "associated 
with" the (say) hypothesis law P^ and an action d A "associated with" the (say) alternative law P A , based on 
the sample path observation X n := {Xi : I G {0, 1, . . • ,n}} of the GWI-generation-sizes up to observation 
horizon n G N. Following the lines of Stummer and Vajda [67] (adapted to our branching process context), 
for BDM let us consider as admissible decision rules S n : tt n h-> {d-u, d A } the ones generated by all path sets 
G n G il„ through 

6n(Xn) ■= fcJXn) ■ = 

as well as loss functions of the form 

L(d%,H) L{d-u,A) \ Li a \ 

L(d A ,H) L(d A ,A) J - \ L H ) [ * U) 

with pregiven constants L A > 0, L-^ > (e.g. arising as bounds from quantities in worst-case scenarios); 
notice that in (80), d% is assumed to be a zero- loss action under T-L and d A a zero-loss action under A. Per 
definition, the Bayes decision rule S Gn min minimizes - over G n ~ the mean decision loss 



dA, 


if X n G 


G n 


d-u, 

■( 


if X n i 

L A 
L H 


G n 

) 



C(S Gn ) := p^ or -Cu-Pr{8 Gn {Xn)=d A 

= Vu Wl ' ^W ' PuA G n) 



u 



p A nm ■C A -Pr(S G „(X n ) = d H 



A 



P P A 0i ■ CA ■ PA,n(^n - G n ) 



(81) 



pP; lor for A. As a 



for given prior probabilities p^' Ior = Pr(H) G]0,1[ for U and p^ nor := Pr(A) = 1 

side remark let us mention that, in a certain sense, the involved model (parameter) uncertainty expressed 

by the "superordinate" Bernoulli-type law Pr = Bin(l,pS 10 ) can also be reinterpreted as a rudimentary 
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static random environment caused e.g. by a random Bernoulli-type external static force. By straightforward 

{prior j °\ 
Z n > p^ior n > leading to the minimal 

mean decision loss, i.e. the Bayes risk, 

Tin := min£((5 G J = £(£<?„ ml J = / min Ip^ 1 " 1 L u,P A im L A Z n \ dP H . n . (82) 

Notice that - by straightforward standard arguments - the alternative decision procedure 
take action d A (resp. d%) if L-h ■ p^ s (X n ) < (resp. >) L A ■ p^ os (X n ) 

with posterior probabilities p^ s (X n ) := pri p " -^^ =: 1 — p A ° s (X n ), leads exactly to the same 

actions as 8o n min - By adapting Lemma 6.5 of Stummer and Vajda [67], one gets for all Lu > 0, L A > 0, 
p^ lm e]0, 1[, A e]0, 1[ and neN the upper bound 

fi n < K\ K l H x H x (P A , n \\P H , n ) , with K n := p% lor L n , A A := (1 - p p £ or )L Al (83) 

as well as the lower bound 

which implies in particular the "direct" lower bound 

. max{l lT ^} . max{l,i^A} 

Tin > -^ H ( x —, ■(Hx(P A AP<H, n )r a * { -^ } ■ (84) 

By using (83) (respectively (84)) together with the exact values and the upper (respectively lower) bounds of 
the Hellinger integrals H\ (P A , n \\pH,n) derived in the preceding sections, we end up with upper (respectively 
lower) bounds of the Bayes risk TZ n . For different types of- mainly parameter estimation (squared-error type 
loss function) concerning - Bayesian analyses based on GW(I) generation size observations, see e.g. Jagers 
[30], Heyde [23], Heyde and Johnstone [24], Johnson et al. [32], Basawa and Rao [4], Basawa and Scott [6], 
Scott [60], Guttorp [19], Yanev and Tsokos [74], Mendoza and Gutierrez- Pena [54], and the references therein. 

Alternatively to the BDM applications above, let us now briefly deal with the corresponding NPT framework 
with randomized tests T n : H n h-> [0, 1] of the hypothesis P-u against the alternative P A , based on the 
GWI-generation-size sample path observations X n := {Xi : I € {0, 1, . . .,n} }. In contrast to (81), (82) a 
Neyman-Pearson test minimizes - over T n - the type II error probability J„ (1 — 7^)dP_4.„ in the class of 
the tests for which the type I error probability J n T n &P-H. n is at most <j £]0, 1[. The corresponding minimal 
type II error probability 

£, (PA,n\\Pu.n) := mf / (l-Tn)dP A , n 

can for all <; s]0, 1[, A <E]0, 1[, n € N be bounded from above by 

MiVHiV) ^ min i (l-A)-(-) • (Hi_A(P^n||i^,n) J , 1 \ (85) 



which is an adaption of a general result of Krafft and Plachky [35], see also Liese and Vajda [47] as well 
as Stummer and Vajda [67]. Hence, by combining (85) with the exact values respectively upper bounds of 
the Hellinger integrals H\-\ (PA,n\\Pu,n) from the preceding sections, we obtain for our context of GWI 
with Poisson offspring and Poisson immigration (including the non-immigration case) some upper bounds of 
£<; (PA,n\\pH,n), which can also be immediately rewritten as lower bounds for the power 1 — £ ? (P A , n \\P-H,n) 
of a most powerful test at level <j. In contrast to such finite-time-horizon results, for the (to our context) 
incompatible setup of GWI with Poisson offspring but nonstochastic immigration of constant value 1, the 
asymptotic rates of decrease as n — >• oo of the unconstrained type II error probabilities as well as the type I 
error probabilites were studied in Linkov and Lunyova [52] by a different approach employing also Hellinger in- 
tegrals. Some other types of GW(I) concerning Neyman-Pearson testing investigations different to ours can be 
found e.g. in Basawa and Scott [5], Feigin [14], Sweeting [68], Basawa and Scott [6], and the references therein. 

For the sake of brevity, a further more detailed discussion of GWI statistical issues along the lines of this 
section as well as power-divergences-connected goodness-of-fit investigations will appear in a forthcoming 
paper. 
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Appendix A: Proofs and auxiliary lemmas 

A . 1 . Tool and proof for Section 3 

Lemma A.l. For all real numbers x,y,z > and all A £]0, 1[ one has 

x x i/- x - (Xxz^ 1 + (1-X)yz x ) < 
with equality iff — — z. 

Proof of Lemma A.l For fixed x := xz x ~ l > 0, y :— yz x > with x ^ y we inspect the func- 
tion g on [0,1] defined by g(X) := x x y 1 ~ x — (Ax + (1 — X)y) which satisfies g(0) = g(l) = 0, <?'(0) = 
ylog(x/y) — (x — y) < y((x/y) — l)—(x — y) = and which is strictly convex. Thus, the assertion follows 
immediately by taking into account the obvious case x = y. □ 

Proof of Formula (30): For the parameter constellation in Subsection 3.3(a5) we employ as upper bound 

for d>\(x), i£Nn the function 

9x[x) • \ 0, ifx>0. 

Notice that this method is rather crude, and gives in the other cases treated in the Subsections 3.3(al) 
to (a4) worse bounds than those derived there. For the calculation of the Hellinger integral, we first set 
e := 1 - e ^(°) e]0, 1[. Hence, we obtain for all n £ N\{1} 

V j cxp{0 A (w„_i } < > j exp{0 A (w„_i)} 

u: n -\— w T1 _i— 

= exp{^A(w„- 2 )} - e = exp{(/5x(wn-2)}- [1 -e-exp{-<^ A (u;„_2)}] 
< exp{^ A K- 2 )-e-e-^^"- 2 )}. 

In the current setup of Subsection 3.3(a5) we have j3 A ^ P%, which means that lim. E _ ) . 00 (/>a(x) — — oo (cf. 
(p-xiii)). But this together with the nonnegativity of (p\ implies swp xeN exp{<f>x(x) — e • e~ Vx ( x '} =: 5 < 
1. Incorporating these considerations as well as the formulae (7) to (12), we get for n = 1 the relation 
H\(PA.n\\pH,n) = exp{0>, (uj n -i )} < 1 and for all n £ N\{1} as a continuation of formula (9) (with the 
obvious shortcut for n = 2) 

oo oo n 

H X (p A , n \\p H , n ) = E---E U z i%) 

cji=0 w n =0fc=l 

oo oo n — 1 

= e- e no*) 

Wi= w n _i=0 fc=l 



i { (f A (oj n ^)) x (/«K-i)) (1 A) - (A/a(w„-i) + (1 - A)/«K_x))} 

e n«(-)- cx p{-/AK-» e [ „ , J ■exp{^K-x)} 



■ cxp ■ 

oo n — 2 



cji— w n _2-0fc-l 0J n —± = 

oo oo n— 2 



^ E-" E II 2 SM- ex P{^(a;„-2)-e-e-^^-)} 

CJi— ixJ„,„ 2 — fc— 1 

oo oo n — 2 

* '-E- E nO«) ^ ^ sln/2i - ^ 

LJi— W n _2— k — 1 

Hence, #a (Pa,„|LP«,„) ™ 0. □ 

Notice that the above proof method of formula (30) does not work for the parameter setup in Subsection 
3.3(a6), because there one has S = sup,j, eNo exp{<p\(x) — e • e~' px ^} = 1. 
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A. 2. Proofs of Section 4 

PROOF of Lemma 4.1 Recall the fundamental nonlinear recursion of (a« ) (cf. (37), (38)), the corre- 

sponding "substitute" inhomogeneous linear recursion of (a n A ) „ (cf. (39), (40), (41)) and its homogenous 
linear relative (ai 9 * ),/l °"\ eNo (cf. (31), (32)) which by (34) and (40) takes the form 



(q*)Mm ,_ (q*),hom _ Jq*),T ( Jq$),h 



[Sl Mm ) =^ t) ^-{d^ T Y), neN, (87) 



with d,( q >- ' ,T = g* • e x o g]0, 1[. By construction, one has 



^ho™ < a {p for ^ n e ^ ag ^j ag Um ^fc*),*™ = Um £*) = Jq*) 



n~ >oo 



t(«x) 



As an auxiliary step, let us compare a; i— > £V A (x) = q* ' eX ~ ft\ with the quadratic function 

T^ a; (x) := ^-e^ ar+q^e^ I 1 - x 9a ; I ■ x + x 9a M 1 - q*e x o + ^~ e x o • 

Clearly, we have the relations Tj'* ) (»J , * ) ) = x^ = fi** ) (»5 , * ) ), ^-(a:?**) = 9*-e^ A> = ^^(x^ 9 *^, 
and — g^ — (x) < — 1^2 — (x) for all x €]x , 0]. Hence, T^ A (•) is on ]x A , 0] a strict lower functional bound 
of £\ (•)■ We are now ready to prove part (a) by induction. For n = 1, we easily see that a-f < a^ 

(9*) r ( X iq *' 1 ) 2 (</*) -i 

iff e x ° • { 2 x A + 1 j — 1 < 0, and the latter is obviously true. To continue, let us assume that 

ai ?A } < ai 9 A } holds. From this, (41), (87) and (88) we obtain 

< $) = fe^tx^.U-e^Y) 2 = ffPLtf^-M?' 



2 



<£? 



< ^ - d^ T ■ alP - xf^ • (l - d^ T ) . 
Thus, a^_L X < a„+i holds. In order to show (b), we make use of the straightforward representation 

alP = £ (d<^)"" 1 "* • UP + x^ • (1 - *«?).*■)) 
which implies that the sequence I a n A ] is strictly decreasing since for all A: € No there holds by (41) 

E ^+ x ^.(i-d^> T ) < tJ*>(o) < #?>(o) = ?*-A < o. 

The final assertion follows immediately from (88) and the closed-form representation (34) with the choices 
K\ K2, x, v, c given just right after (42). □ 



Proof of Lemma 4.4 For "Psp\('Psp,3d U Vsp,a) we deal with the fundamental nonlinear recursion of 
(on ) , G £ {E,U} (cf. (48), (27)), the corresponding "substitute" inhomogeneous linear recursion of 

(tin ) (cf. (49), (50), (51)) and its homogenous linear counterpart {tin ° m )„ FN ( c f- (31), (32)) which 

by (34) and (50) takes the form 

tif )Mm :=Q, a^ )Mm := ^• S (aS Mm ) = a# ?) (l- (d^*)") , n € N, (89) 
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with d(9?)< s = 1 - 9x ~a } x G ]d^' T , 1[. By construction, we obatin 

_(,?),/.<m. = fl (g«) ) -(9?),^™ > a ^) for ^ n G N ^ 1}; and Hm -foghorn = ^ a (g?) = s (^) _ (go) 



n— >-oo 



In analogy to the Proof of Lemma 4.1, we use the quadratic function 



V X ° 7 — Ca°) 

which satisfies T^O^) = x^ ] = ^\x^ ] ), ff\o) = qf - /3 X = £ a 9a) (0), and ^Jjf-fc) < 

— 1^3 — (x) for all x GJXq A , 0]. Hence, T A A (•) is on ]xq , 0] a strict upper functional bound of £ A A (•). To 
start with the proof of part (a) , let us first observe for n = 1 the obvious relation af x = q x — P X — a^ =0. 
Furthermore, let us assume that a„ A > a„ A ' (n £ N) holds. From this, (51), (89), (90) and the appropriately 
adapted version of a n ° m we obtain the desired inequality a£±.l > a„+i by estimating 



(4 9?) )' 



> Ci 9?) («i ??) ) - d^> S ■ af x) - (qf - &) > a\H - d^- S ■ a^ (qf /3 A ) . 



Moreover, the property (b) follows from the representation 



n — l . , 

^ = £ («**>■*) -(74^ + (a? -a)) 

fc=0 

which implies that the sequence ( ol? A ) is strictly decreasing since for all k € No one has /5/f* + (q x — 

V / riGN 

Px) < 0. Finally, part (c) follows immediately from (90) and the closed-form representation (34) with the 
choices K\ K2, x, v, c given just right after (52). □ 

A. 3. Proofs of Section 5 

Proof of Theorem 5.1 As already mentioned above, one can adapt the proof of Theorem 9.1.3 in Ethier- 
Kurtz [13] who deal with drift-parameters i] = 0, k, — 0, and the different setup of a— independent time-scale 
and a sequence of critical Galton- Watson processes without immigration with general offspring distribution. 
For the sake of brevity, we basically outline here only the main differences to their proof; for similar limit 
investigations involving offspring/immigration distributions and parametrizations which are incompatble to 
ours, see e.g. Sriram [61]. 

As a first step, let us define the generator 

2 
A.f(x) := (r? - k. • x) f(x) + y • x ■ f"(x), f e C c °°([0,oo)) , 



which corresponds to the diffusion process X governed by (61). In connection with (58), we study 
Ti m) f{x) := EP. 



( 1 / rax 



( m ) I y( m ) 



, x e £H : = In 0j /eC c °°([0,oo), 
m 



where the Y Q \™ , Y" are independent and (Poisson-/3. respectively Poisson-a, ) distributed as the mem- 
bers of the collection y(' n ) respectively Y^ m \ By the Theorems 8.2.1 and 1.6.5 as well as Corollary 4.8.9 of 
[13] it is sufficient to show 

lim sup o- 2 m(T<r ) f(x)-f(x))-A.f(x) = 0, / e C c °°([0,oo)) . (91) 

But (91) follows mainly from the next 
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:>,() 



Lemma A. 2. Let 

/ n 

with the usual convention Sq := 0. Then for all m E N, x G E^ m > and all f £ C£°([0, oo)) 



,(m) _ a(rn)\ ^ y(m) _ Q (m) \ „ g Nj m g JJ, 



e lmJ (a:) := PP. 



(Sf) -(1 - v) f" ?i m) x + ^-+ vJ^StJ /"(*) dv 



o 2 m-(T { . m) f(x)-f{x))-A.f{x)\ + R (m \ where lim P (m) = 0. 



(92) 



Proof of Lemma A. 2 Let us fix / e C£°([0, oo)). From the involved Poissonian expectations it is easy 
to see that 



lim 

m— »oo 



a 2 m(T. (m) /(0)-/(0))-A./(0) 



, 



and thus (92) holds for x = 0. Accordingly, we next consider the case x € P( m '\{0}, with fixed mgN. From 



PP. 



5, 



(m) 



= 0("») + o£L weobtain 



pp. 



(S^) 2 xf"{x)f (l-v)d 



= \[p { ™ ) -x+^)f"{x) =: a 



Furthermore, with b mx := b := a + ^fxfm ■ S%?J = ± (X)fefi *ofc + *0 ) we § et on { S rnx 7^ 0} 



(m) 

■ J (X) =: a mx — — - =: a — 
/ 2 2 

( m ) ,v-( m )^ ™ o- Q + ™ ro( m ) 



™ - a m. (93) 



1 / rv(™) 

II R {m)„ a 



/" / 9W x + i5_ +t>1 /* 5 Cm) * 



771 1 



x q(" 1 ) 

*~>m a; 



/"(i/)dy 



m f'{b) - f'(a) 

c (m) 



(94) 



as well as 



\ m V m / 

m /'(6) m f(a)-f(b) 






yf"(y)dy-a f f"(y)dy 



X q( m ) x 



foWV 

1 Omx J 



(95) 



With our choice /3. =1 ^- and a, = /3. • -5f , a Taylor expansion of / at x gives 

f(a) = f(x) + -L.f(x)(pi m) -r 1 -K.-x) + o(-), 

am \ / \ m / 

where for the case r\ = n = we use the convention o ( — ) = 0. Combining (93) to (96) and the centering 



(96) 



PP. 



S 



(m) 



0, the left hand side of equation (92) becomes 

S m ) 



EP. 



(S&0) x(l - v) f" pi m) x + ^- + v J?- S %> - f"(x) dv 



EP. 

1 



^■S£J-(f(b)-f'(a) 

„ (™) \ 

,(m) 



EP. 



^■StJ-f'(b) + m-(f(a)-f(b)) 



Z \ m 



1 / (m) \ 

(ep. [m] -/(«)) - 5(' 3 -""^ + ^rj^» 



m ■ < PP. 

1 



a=i 



/(*) - -^/(s) 



(?7 - k. ■ x) - P. ■ r) + n. ■ x ■ f'(x) + - 



1 _p(m)_a. 



(m) 



fix) -m-o[ — 
m 
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which immediately leads to the right hand side of (92). 



31 

□ 



To proceed with the proof of Theorem 5.1, we obtain for m > 2k. jo 1 the inequality /3. > 1/2 and 
accordingly for all v <G]0, 1[, x <G E^ 

(m) i (m) ( rax \ - 

/3^+^+^^ ) = (l--)-^4 m) + (l--)^+^E F oI ) +?oJ > *• V- 

Suppose that the support of / is contained in the interval [0,c]. Correspondingly, for v < 1 — 2c/x the 
integrand in e' - m '(x) is zero and hence with (97) we can estimate 

(t») 2 x(l - *) (V (V" } * + ^ + vM 5W] - /"(*)) d« 



< 



fe) 2 ^(i-^)-2ii/"iL^ < ^•te) 2 fiA-) 2 |i.rii c 

V(l-2c/a0 V ' K ' \ x J 



From this, one can deduce lim m _>. 00 sup,j, e £( m ) e^ m '(x) = - and thus (91) - in the same manner as at the 
end of the proof of Theorem 9.1.3 in [13] (by means of the dominated convergence theorem). □ 



The following lemma is the main tool for the proof of Theorem 5.3 below. 

Lemma A. 3. Let (kj., K-h, n, A) £ {Vni U7 ? sp,i)x]0, 1[. By using the quantities k\ := Xka + (1 — ^) K u > 
and A\ := ^/Xk\ + (1 — X)kiL > k\ from (64), one gets for all t > 

(a) lim m-(l-^ Tn) ) = ^ > 0. 

(6) Jimm 3 • fl « = _ Hl-X)( K j- K uf = _*t_£ < Q . Um m . (1 _ /3 M ) = ka > a 



(c) lim m ■ Xq 



(m) 



2a 4 

Aa - «A 



< 



2a 4 

,• 2 r (m) (^ - «a)' 

lim m ■ l y ' ~ — : 

Tl— >00 2(7 



> 0. 



(d) lim m-(l-d(™)' s ) = ^^ > 0. 

(e) lim m • (1 - d (m) ' T ) = ^ > ; lim m 2 . x { ™ ] ■ (1 - d (m) ' T ) 

771— J-OO (J ??l— J-CXD 

(/) lim m-(l-d^' S d^ T ) = 3Aa+ / A > 0. 

777,— >OG 2(7 



A A • (A A - k\) 



< 0. 



( 3 ) lim (d {m} ' s ) 

771— J-CXD V / 



Aa + ka 
cxp <i - • t )■ 



/ \ a rat 

(h) lim (d {m) ' T ) = cxp{-AA-t} < 1. 

(i) lim U^d^Y" = exp(- 3AA + ^ 



n < i 



iy^mt\ 2ct 2 • A 



A 



V ( Aa - «a 
o L k 2 ™*J ' 4 CT 2 



WJm^&Ll ^ T^-(^r^) "(1-e-^r > 0. 



A 



-A 



m lim m 7 { " l) - ( A ^ ~ K *)" 
(l) lim m-C,, 2 mt i - = 

m— >oo L J cj^ 



-^(Aa+ka)-* _ e -Av* e -|(Ax+Kx)-« (l _ e ~ A A-t~) 



(m) lim z? L ™2 mtJ 



V (Aa - K X y 
a 2 Aa 



A A - k a 2 • A A 

1 — e -3(3Ax+Kx)-t e - A A'* — e -j;(AA+K A )-t 



> 0. 



3A A + k a 



Aa - k a 



> 0. 



Proof of Lemma A. 3 For each of the assertions (a) to (m), we will make use of l'Hospital's rule. To 
begin with, we obtain for arbitrary /j,, v £ M 



lim m 

m—toc 

lim m 

m— >-oo 
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i-(/?S m) n/4 m) r 



n ■ {^y-\^yJ!A + v ■ (^^(^j-i-pL 






K-H 



From this, (a) follows immediately and (b) can be deduced by 



1- 2 (m) ,. m 

lim m ■ a\ = lim — - 

n—toc to— >-oo 2(7 



A-M^-cejrV- 1 ^) 1 ^) 



(to)nA/o(to)n-A 



+ (1-A).«k1-(^)W) 



A(l-A)(^-^) s 
2a 4 



•T2 



A m ) 



A m ) 



A™) 



„( m ) 



For the proof of the first part of (c) , we rely on the inequalities Xg < x < x Q (m S N) , where x 
and Xq are the obvious notational adaptions of (46) and (56), respectively. By using (a) and (b), one can 
calculate 



lim m • Xq = 

m— >oo m— >oo 



lim ( q ( ™' 



/1 V rrL )\ 

m - (1-9A )- 



(to) 



)-y^To 



(m). 



(m) 



i-?r) -2. 9 r-^ 2 -«i 



Aa - «A 



,M 



From (46), (a), (b) and linx^^oo pi * =1 we obtain the limit 



lim h 



(<T) 



= lim max < — P\ 



A m ) 



1 (to) 



= lim 



2 (to) 

m ■ a) ' 



(m) 

™- [l-q\ 



which implies 



lim r 

771— >00 



("») 



lim ■ 

m— ?oo 



9a 



m • (1 - g A 



(to)^ 



/Y fi (m)x\ 2 o ft (9i m) ) (to) 9 (to) 



= 0, 



Aa - k\ 



and thus the first part of (c). The second part is an immediate consequence thereof. Assertion (d) follows 



from (b) and (c) by 



lim m • (1 - d (rn) ' s ) 



A m ) 



lim 



m ■ x 



(to) 



Aa + k\ 
2a 2 ' 



For the first part of (e), we use the general limit lim x ^o - — = 1, to get with (a) and (c) 



lim m ■ (1 - d (m )' T ) 



lim m • ( 1 

771— >OG \ 



9a 



9a 



A m ) 



e o — l 



„(to) 



Aa 

^■2 



From this and (c), the second part of (e) is obvious. The limit (f) can be obtained from (d) and (e). The 
assertions (g) respectively (h) respectively (i) follow from (d) respectively (e) respectively (f) by using the 
general relation lim m ^oo (l + ^-) = e hmm ^° CiXm . The last four limits (j) to (m) are straightforward impli- 
cations of (a) to (i). □ 



Proof of Theorem 5.3 It suffices to compute the limits of the bounds given in Corollary 5.2 as m tends 
to infinity. This is done by applying Lemma A. 3 which provides corresponding limits of various involved 
quantities. Accordingly, for all t > the lower bound (66) can be obtained from (62) by 



lim exp < Xq 
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x, 



(m) ?7 



ff 2 ]_ _ d (m), 



1- fd< 



33 



(m) V 



L ° O 2 



^•l- 2 h + cx tl -xr + c= J 



la 2 mt] 



j(m) 

-Lcr 2 mtJ 



,. I (m) 

= lim exp < to • ccq 



X, 



(„,) 



v d^ T 



a 2 m ■ (1 - d(™'< T ) 



exp ■ 



+ 771 ■ X 
Aa - «A 



(m) V [° 2 rnt\ 



X, 



a 2 



rj a 



-Y<j-mt\ m 



l- s 



(m) 



[a 2 mt\ 



a( m ) 
^Lff 2 mtJ 



a 1 A A 



(Aa - K X f . t , 
2a 2 . Aa li 



fl _ e -Axt) _ ^Z^ . ^ . ^ 



-A A -i\ 



Xn 



er^ er^ 



7] /Aa-k a 



4a 2 



Aa 



1-e 



-Aa-M 



exp 



Aa - ka 



x °-t 



(1 - e-*.'<) - 4j (Ax - «») • t + 4'V) • Xc + ^ • Lf(() 



For all £ > 0, the upper bound (66) follows analogously from (63) by 



lim exp < Xq 



(m) 



y( m ) ]_ 

o a 2 ' 1- dM< s 



(l - (d^ s ) 



\a 2 mt\ 



+4 m ^-w - c ( L ;Lj-^r 



(m)_?7_ 
<7 2 



(m) S (™) 

[cr 2 mtj 



lim exp < m ■ x 



771 cr 2 



d(m),S 



to- (l-dH.s) 



fl - (V m >' s ) 



[<r 2 mtj 



+ 771 • CCA -^ 



(m) V [^ 2 mt\ 



(m) AT 



(m) 



/?? 



m ' CL<r 2 mtJ ' 



m 



j(m) 
'Lo- 2 mtJ 



exp ■ 



Aa - ka 



exp ■ 



An 



e -±(A A + KA )*^ Aa ka r ? 



(Aa - k\Y 
a 2 

n (A A - k\Y 
<J 2 Aa 

Aa - k\ 



7] 8cr 2 
cr 2 ' A A + k x 

e -5(AA+K A )-t _ e -Ax-t e -5(A A + KA )-i A _ e -A A -n 

A A - k x 2 • A A 

1 — e -ij(3A A +K A )-t e _AA '* — e -5( A A+fA)-t 



(7 Z t 



X 



3A A + k x 



Aa - K X 



Xq j- 



'/ 



(A A + k x ) 



(l_ e -|(Ax + ^_^ (AA _„ A) . f 



U^^-Xo - ^-C/f(t) 



n 



^4.^. Proofs of Section 6 

We start with two lemmas which will be useful for the proof of Theorem 6.1, and which can be easily seen 
by induction. They deal with the sequence (a„ ) from (17). 

Lemma A. 4. For arbitrarily fixed parameter constellation {Pa, PH-i a A-i a H) G "P , suppose that qx > 
(A s]0, 1[) and liniA/'i <?a = Pa holds. Then one gets the limit 
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VneN: lim ai qx) = 0. 



34 

(97) 



Lemma A. 5. In addition to the assumptions of Lemma A. 4-, suppose that A i— > q\ is continuously differen- 
tiable on ]0, 1[ and that the limit I :— lim,\/<i -Sr- is finite. Then one gets the limit 



VneN 



.. dci 9x) 
urn — — — 
\yi dX 



l+Pu-P; 



= u n := 



(1-(Aa)"), ifP A ^l, 
if Pa = 1, 



n-(l + p n -l), 

which is the unique solution of the linear recursion equation 

u n = I + Ph-Pa + PA-Un-i, u = . 

Furthermore, 



VneN 



Elim 
X/-1 



= E u ^ = 



l + fi-H-tiA 

i-Pa 



k=l 



A/'l d\ 



fe=l 



i-Pa 

^■{l + fiu-l) 



(i - (PaT) 



if Pa + 1, 
if Pa = 1. 



We are now ready to give the 

Proof of Theorem 6.1 

(a) Recall that for the setup (Pa, Ph, a A, a u) € ("Pni U Psp,i) we chose the intercept as p>, := pf := a^a^ 



a „,i-A 



and the slope as q\ 
Ea 

'/A 



?A 



E ._ 



PaPh i w hich in (20) lead to the exact value V\_ n of the Hellinger integral. 
Because of —P\ — ct\ = as well as liniAyi q\ — Pa, we obtain by using (19) and Lemma A. 4 for all n € N 



limF A ,„ := limexp ol to » ■ uj + Y b^ x ' qx) } = lim exp { a<?*> ■ u + -± V u , 



(9X) 



A/ 1 



which leads by (71) to 



I(PA,n\\P-H,n) = hm 

A/ 1 ! 



1--HA(-P>l,n||i^,n) ,. 1-Vx,„ 

= lim 



A-(l-A) 



a/<i A- (1- A) 



= lim 



-Vx, n 9 



A/-1 1 - 2A <9A 



(<7a) i P X 



yA fc=i 



feO 



LngN, 



= lim 

A/-1 



da 



(«a) 



(9A 



w 



j9_Pa 
<9Ag A 



£4' 
fc=i 



<Za A? dX 



fc=i 



(98) 



For further analysis, we use the obvious derivatives 
dp\ 



dX 



. a A \ d p\ p\ (a A Pn\ dq\ ( p A 

PA log , jrr— = — log — , — - = q X log — 

a H J dXqx qx \ctuPA dX \Pu 



(99) 



where the subcase (Pa, Ph, o> a , au) G Vta (with p\ = 0) is consistently covered. From (99) and Lemma A. 5 
we deduce 



lim 



da 



(<3a) 



Vi <9A 
and by means of (97) 



•w 



(Pa log (f£) - (Pa - Ph)) ■ '-^f- ■ ^o if Pa ± i, 
n-(^log(|^)-C8^-i8«))-a;o if#A = l, 



VneN: lim 

VI 



5 M f n fc) 

d\qxl^ ak 



fe=i 



= 0. 



For the last expression in (98) we again apply Lemma A. 5 to end up with 







lim — • > — -i 
A/-1 ?A ^ <9A 



,(<7a) 



<^a\Pa^&{^)-(Pa-P-h)] 
Pa{i-Pa) 



0, 



1~Pa 



(1 - (P A T) if fo ? 1, 



n.(n + l)^-- ^log(|^)-(^-/3 H ) 



iffti = l, 



(100) 



which finishes the proof of part (a). To show part (b), for the corresponding setup (Pa, P~H, a A, a u) G 
'Psp\'Psp,i let us first choose - according to the Section 3.2 - the intercept as p\ := p\ := a A a^ and the 
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slope as q\ := q^ := P^Pn , which in part (b) of Proposition 3.4 lead to the lower bounds B^ n of the 
Hellinger integral. This is formally the same choice as in part (a) satisfying livax/^ipx = a a, liniA/i q\ = Pa 
but in contrast to (a) we now have — px — a\ ^ but nevertheless 



q a = 0. 
From this. (19), part (b) of Proposition 3.4 and Lemma A. 4 we obtain 



lim — j3 x 

vi qx 



lim Bf _ 

xyi A - n 



lim cxp < a;; 

xyi 



(qx) 



Wo 



Px_ 

qx 



fe=l 



(?a) 

fe 



PA 
?A 



j3x - ax 



(101) 



and hence 



I{PA,n\\ P H,n) < Jim 



1 



dL p?L 

= lim 



S 



x/iA'(l-A) A/<il-2A <9A 



a^^o 



lim 

A/>1 



<9A 



-wo ■ 



_9_PA 

dXqx 



E< 

fe=i 



(9X) , PA 

fc "T" „ 

<7A 



E 






^ 



8 



OX 



dX \q\ 



PX n 

— Px - ax 

qx 

Px 



Px - a\ 



(102) 



In the current setup, the first three expressions in (102) can be evaluated in exactly the same way as in (99) 
to (100), and for the last expression one has the limit 



d fpx a 
aV — Pa - "a 
oX \qx 



V.i 



Px , 
— loj 

qx 



faApH 

\anPA 
(oaPk 
\anPA/ 



Px + (Pa- Pn) - [aA - ay) 

qx 



log (^L\ 



Pn 
Pa 



ay 



which finishes the proof of part (b) . 



□ 



Proof of Theorem 6.3 Let us fix (Pa, Pn, a a, an) S PspXPspa, wo € N, n £ N and y e [0,oo[. The 
lower bound E^^ an of the relative entropy is derived by using as a linear upper bound </>^ (cf. (15)) for (px 
(A €]0, 1[) the tangent line of <j>\ at y. This corresponds to <fix(x) := (p^ — ax) + (q^ — Px) x (x £ [0, oo[) with 
Px ■= Px(y) ■= 4>x{y) - y<t>'\(y) + ax and qx := q\(y) := (f>' x {y) + Px, implying q x > because of (p-xii). As a 
side remark, notice that this <fix(x) may become negative for some x £ [0, oo[ (which is not always consistent 

with goal (Gc) for fixed A, but leads to a tractable limit bound as A tends to 1). Analogously to (101) and 

V 

A,: 



(102), we obtain from (19) and (21) the convergence Ymvx/^i Bx„ = 1 and thus 



I{PaA\Ph,u) > lim 



dX 



kfttV+Zt' 



ox 



9 [PXr, 

oX \qx 



fc=i * fe=i 

As before, we compute the involved derivatives. From (10) to (12) as well as (p-xii) we get 



(103) 



dpx 
dX 



VJ 



( My) 
\fn(y) 

+ Puy 
a^tlog 



/«(y)iog 



My) 



Pav 



( My) 
\Mv) 



A-l 



My)\ x n wr ( Mv) \ x , , 
My) ) - {1 - X)PHy \My)) l ° g 

( My) \ y ■ (aAPn -an Pa) 

\Mv)) My) 



and 



dqx 
dX 



"5T = Pa 



fA(y) 



A-l 



.My) 
+ (1 - X)p n 

Pa (l + log 
Combining these two limits we get 



xp A 



A-l 



A/7 



( fA(y) \ 

\My)J 

My) \ 

My)) 



( My) 

\My) 



*m 



lot 



- XPaV 

My) \ 

My) J 



( My) \ 
\My)J 



My) 
My) 



A-l 



log 



My) 
My) 



(104) 



-Pn 



My) 
My) 



Pn 



My) 
My) 



=: I. 



(105) 
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:',(, 



d_ 
ox 



Px 

qx 



Px - ax 



<1\ 



dX 



PX 



dX 



A/-1 



(qx) 2 
y ■ (oiaPu - cx-hPa) 



px_ 
qx 



fn(y) 



-oca 1 



(Pa - Pn) - (<*A - <*h) 

olaPh 



Pu1a(v) \ 

PaMv)) 



Otn 



Pa 



(106) 



The above calculation also implies that lim,\yi I ^-— ) is finite and thus lirnxyii I Jj — ) X)fc=i a k = by 

means of Lemma A. 4. The proof of I(PA,n\\pH,n) > Ey'^ n is finished by using Lemma A. 5 with I defined 
in (105) and by plugging the limits (104) to (106) into (103). 

To derive the lower bound E k ^ ec (cf. (76)) for fixed k £ No, we use as a linear upper bound </>^ for 4>\(-) 
(A £]0, 1[) the secant line of <j>\ through the points k and k + 1, corresponding to the choices p\ := px(k) := 
(k + 1) • 4>x(k) - k ■ 4>\(k + I) + a x and g A := qx(k) := 0a (fc + 1) - 4>x(k) + /3a, implying q\ > because of 
(p-xiii) and (p-iv). As a side remark, notice that this (j>^(x) may become negative for some x £ [0, oo[ (which 
is not always consistent with goal (Gc) for fixed A, but leads to a tractable limit bound as A tends to 1). 
Analogously to (101) and (102) we get again liniAyi B Xn = 1, which leads to the lower bound given in (103) 
with appropriately plugged-in quantities. As in the above proof of the lower bound Ey-^ an , the inequality 
-f(-P4,n||-P"H,n) — ^fc iT^ follows straightforwardly from Lemma A. 4, Lemma A. 5 and the three limits 

dpx ffA(k)\ X ,, wiu . ni (fA{k)\ f f A (k+l) \\ f fA{k+l) 

ffA(k+iy 



A/-1 



Mkmi)k ,(L&L) 



f A (k+l)klog 



\fn(k+l) 



dqx f fA(k+l) \\ . ffAk+m ffA(k)\ 



OX 



*<*"«(*§) 



^^Mm)-^»<m) 



=: I 



and 



Px 



d_ 

dX \q> 
xyi 



Px - a x 



qx 



ox 



Px 



t)X 



Px 



^HW>) 



k+l 



(qx) 2 

OA 

Pa 



Px H {Pa- Ph) - ("a - a H ) 



qx 
jU(fc+l)log 



( fAJk- 



\fn(k+l) 



OLA 

Pa 



uaPh 
Pa 



an- 



To construct the third lower bound E^ hor (cf. (77)), we start by using for each fixed A s]0, 1[ as an upper 
bound of (j)\ the horizontal line through the intercept sup xGNo 4>x(x). For 'Psp.Sab U 'Psp.Sc, this supremum 
is achived at the finite integer point z* x :— argmax 2 gN 4>x{x) (since liuix^ao 4>x(x) = — oo) and one has 
4>x(zx) < which leads with the setup qx = Px, Px = 4>\(z\) + °^ to the Hellinger integral upper bound 

1 — B u 

B Xn = exp {<p\{z x r ) ■ n} < 1 (cf. (47)). To compute from this the required lim,\/<i A / 1 _ A A ? is not straightforward 
since in general it seems to be intractable to express z* x explicitly in terms of A. However, since lim^/'i <px(x) — 
for all x £ [0,oo [, we obtain by l'Hospital's rule 



lim 



Mx) 



V"i 1 - A 
Accordingly, let us define z* := 
the maximum exists since lim T 



(oca + Pax) 



-log 



a a + Pax 



argmax xeNo I {a a 
,ola + Pax] 



{(' 



+ Pax) 
log 



1 



Pux 



- (a H + P H x). 



an+Pnx 



1 



- (an 



- (an 
Pnx)} = 



Ph x ) \ (note that 
— oo). Due to con- 



tinuity of the function (A, x) H> /_ A , there exists an e > such that for all A e]l 
Applying these considerations, we get with l'Hospital's rule 

1 - exp{0A(^*) ■ n} 



e, 1[ it holds z x 



I(PaA\P' 



.„) > lim 



A(l-A) 



fA(z*) 



log 



Ia(z*) 

fu(z*) 



-1 



+ Iu(z*) 



n > 0. (107) 



In fact, in the current parameter constellation 'Psp.Sab U 'Psp.Sc we have <j)x(x) < for all A g]0, 1[ and all 
x e N which implies Ja(z*) ^ fn(z*) by Lemma A.l; thus, we even get E^' hor > for all n £ N by virtue 

of the inequality 



log \fA(z*) 



> 



f-HJz*) , 1 

7aW) +i ' 
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For the case"Psp,2, the abovementioned procedure leads to z^ = = z* (A g]0, 1[) which implies (j>\ (z^) = 0, 



1-Bk 



B^ n = 1 and thus the trivial lower bound E^' hor = limA/u x(i-x] = follows. In contrast, for the case Psp^d 



one gets z 



A ~ P-h-Pa 



z* (A g]0, 1[) which nevertheless also implies (j>\(z%) = and hence E, 



L.hor 



0. On 



fsp,4j we have sup a;gNo <fi\(x) = <fi\(oo) = and hence we set E^' hor := 0. 

To show the strict positivity E^ > in the parameter case "Psp^, we inspect the bound E '^ ec . With the 
help of a := a, := a. a = a H (the bullet will be omitted in this proof) and the auxiliary variable x := ip- > 0, 



E, 



L,sec 
0,n 



I-Pa 



the definition (76) respectively its special case (79) rewrites as 

.{a + f3 A )-log(?£^f)+P A {x-l) 
~JA^I) (-(« + M ■ log (*ftf j + Pa{* - 1)) 



w 



1-/3.4 



T-iL.sec/ \ 

^0,n (*) = S 



.E- 1 



(a + 1) • log (^ 

'( a + i) . log (0) - x + 



+ (wo + f ) • n] 



if/3^1, 



X0A 



(108) 



To prove that E^ ec > for all w„efj and all n G N it suffices to show that E%£ ec (l) = (|£ H (1) = 



and 



rE, 



L,sec 



we obtain 



(x) > for all x s]0,oo[\{l}. The assertion E ^ ec (l) = is trivial from (108). Moreover, 



Tx E ^ ]{x) 



q+Ps 



a+p A x 



a- 1- 






1-(M" 
/3.A 



W 



1-0.4. 



1-/3/ 



• n, 



1 _ "+1 
a+a; 



Li • *' + {"o - f ) ■ n\ , 
which immediately yields I ^i5 ^ ec ) (1) = 0. For the second derivative we get 



if/3.4^1, 
11/3.4 = 1, 



(109) 



<9 2 



(c*+Pa)-P a i-(M r 

(a+[5 A x) 2 ' 1-/3,4 



W 



i- Pa 



a+fJA 



P't 



■ (a+/3 A x) 2 1-/3.4 



■ n > 0, 



if 0^1, 



(110) 



q+l 



[f •n 2 + (w -f)-n] >0, if/3^ = L 



where the strict positivity in the case /3^ 7^ 1 follows immediately by replacing wo with 1 and by using the 
obvious relation 



i-Pa 



1-/3" 



i-Pa 



drr EVo (1 - #) > 0. 



For the constellation Vspa with parameters /3 



Pa = Pui &A 7^ Q^, the strict positivity of 



E^ > follows by showing that Ey'^ an converges from above to zero as y tends to infinity. In fact, there 
holds lim^^oo y ■ Ey£ an e]0,oo[. To see this, let us first observe that by l'Hospital's rule we get 



v , (aA + Py 
hm y ■ log — — 



/3 



a A + Py 
as well as lim y • I 1 

y^cx) a-u + py 







From this and (75), we obtain lim^oo y ■ Ey£ 



( aA -a n ) 2 . n > in both cases /3 ^ 1 and /3 = 1. 

7,L,ta 



Finally in the parameter case Psp,3d we consider the bound Eyl ° n , with y* — a „ A _" K . Since a^ + /3_4y* = 
a-H + Puy* holds, it is easy to see that El °™ = for all neN. However, the condition I -^-Ey'^ an ) (y*) 7^ 
implies that sup y>0 Ey^ an > 0. The explicit form (78) of this condition follows from 



^ E vi an 



(a A p n -a n P A r 1~(Pa) 

f A (y)(fn(y)V i-Pa 

_l_ o-aP-h — Q'-hPa 
"*" Uu{y)Y 

(u A Pn~an) 2 

U(y)(fn(y)r 



Wo 



i-Pa 



P A (i-P A )fA(.y) 



a A Pu-aHP A 
P A 



if/3^^1, 



\ Qi A 

L 2 



(<* + ¥)•»] - ^ZTrJ' -". if ^ = 1 - 
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□ 



The next lemma (and parts of its proof) will be useful for the verification of Theorem 6.6: 

Lemma A. 6. Recall the bounds on the Hellinger integral m— limit given in (66) and (66) of Theorem 5.3, 



-«/ 



r(i) 



in terms of L\ (t) and U^ it) (i = l,2j defined by (67) to (70). Correspondingly, one gets the following 
A— limits for all t £ [0, oof: 

(a) for all ka g]0, oo[ and all k% € [0, oof with k^ =/= k-h 



\/-\ dX 



(b) for ka = and all k-u E]0,oo[ 



(2), 



r(i)/ 



hm — A --^- = hm 



A/1 



9A 



a/i dA 



hm |-A2 



lim 



lim 

A/1 



dX 

dLf{t) 



2a 2 ' 

.2 4-1 






<9A 4 



.. du^>(t) ,. ac/r ; W 

hm £— ^ = hm 



A/i 



dX 



A/1 



9A 



0. 



0. 



(Ill) 

(112) 
(113) 
(114) 



Proof of Lemma A. 6 For all k^, kh € [0, oof with k^ ^ k-h one can deduce from (64) as well as (67) to 
(70) the following derivatives: 



dL{ 1] (t) 
dX 



dX 



1 J t / Aa-Ava 
2<t 2 1 2 I A A 



{ K A ~ 



?t)[2e 



-2Axt 



,-Aii 



-A A t 



1 



-Ax* 



Aa 

lJ Ax-nx 
41 A A 



' K .4 ~~ K « ~~ 2A a (k^ - Kfi)) 



Ax 



Ax - K\ 



Aa 



2 „2 



' H 



1-e 



-Aa* 



A 



A 



K A~ K u~ 2A a (k^ - k w ) - 



Aa - « A 
A A 



_ K w) 



,(115) 



«9A 



^.fA^^y.i-^.^ 



1_ I 1 A A -« A 
rr 2 1 2 ' A A 






A 



Af 



v Mj 



4) "5 e-^ (AA+KA)t (4 - 4 + 2A A (^ - «„)) 



KA\+kx)< . 



-Axt 



Aa 



{ k a ~ K u + 2A A («u - Ku)) 



I e -i(Ax+re A )t _ g-A\t 



A A 



«3l _ K W - 2A A (Kyl - K W )) 



+ 



l (A 



K-X 



Ax 



e -|(A A+KA )t (^ _ ^ + 2Aa( ^ _ ^ 



t e _ i (3AA+KA)t (3 ^ _ ^ + 2Aa( ^ _ Kk)) + e -i(Ax+«A)t 



1-e" 

AT 



-Ait 



"«J 



dUi 2 \t) 
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OX 



(AA - Kxf 

A a (3A a + k a ) 



t 



= -i(3A A +K A )t 



.K~ 



K 



H 



2A 



ka - mi 



-X 



e 2 



i(3A A +K A )i 



A A - K X 
Ax 



e 2 



3A A + k a 



(A A +K A )t 



V H 



2A A 



KA - K-H 



K 



n 



e 2 



|(A A +K A )t 



-A x t / 2 



A 



-A 



K 



2A A 

2 



K.A - K H 



te 



-A x t ^A 



K A Kj 



H 



2A, 



H 



2A 



A 



HA + K n 



' H 



1 

A7 



2A A 

A A - nx 
3A A + k a 



KA + K H 



Ax K\ Kj^ — K-^ 



A| 



l-e - 3 



i(3A A +K A )t 



") _ e -^(A A + « A )4 + e -A A t 



(116) 



If k^ g]0, oof and k^ S [0, oof with kj± =^ re^, then one gets lim A /i 1 A A = limwi k a = k_4 > which implies 
(111) from (115) to (116). For the proof of part (b), let us correspondingly assume kj, = and ku g]0,oo[, 
which by (64) leads to k a = Ku ■ (1 — A), A A = k% • \/l — A and the convergences lim A yi A A = \ivi\^i k\ = 0. 
From this, the assertions (112), (113), (114) follow in a straightforward manner from (115), (116), (116) - 
respectively - by using (parts of) the obvious relations 



lim 

A/U 



lim — - = 0. 
Vi A A 

l_ e -cvt 



lim 



A a ±k a 



t 



Ax 
for all c x E {A x 



lim 



A A - k x 



1 



Vi A A + k a 

A x + kx 3 A A + k a 



c A L 2 2 

In order to get the last assertion (114) we make use of the following limits 

13,. 4k-h 



lim 



= lim 



A/-1 A A - kx 3A A + kx xyi (k u - K U ■ \/l- A) • (3 k h + k h ■ \/l -A) 



3k-h 



and 



lim -— 
Vi A A 



I _ e -^(3A A +K A )i ^ _ e -A A t ^ — e -5(A A + K A )i 



3A A + K X A\ - K X Ax - K X 

To see (120), let us first observe that the involved limit can be rewritten as 



= 0. 



1 



a™ ] A a (A a -k a ) 



1 1 



e -i(3A A +K A )t + g -A A t _ e -I(A A + KA )i 



1 - e -^(3A A +K A )t 

A^ 



1 



3A A + kx 3(A A - kx) 



(117) 
(118) 

(119) 
(120) 

(121) 
(122) 



Substituting x := yl — A and applying l'Hospital's rule twice, we get for the first limit (121) 



1 _ 1 
lim - 

x\,0 



e -^-(3x+X 2 ) 1 e -K-HtX _ e ~^-(x+X 2 ) 



K H ■ (x 2 — X 3 ) 



lim 

x\,0 



^(3 + 2s) e-^ 1 ^^ 2 ) - K H te- K ^ tx + ^ (1 + 2x) e -^ i ( x+x2 ) 



rb-X/ " I ^3j OiX I 



lim - 

x\0 

1 



^(3 + 2x) 2 + 



3 



'-(3x+x 2 ) + K 2 f 2 e 



2 j.2 e3 — K,-i-itx 



n ' 



£f(l + 2xf -K U t 



^-{x+x 2 ) 



k 2 , ■ (2 — 6x) 



tjrXinj L 



3 



«- 2 / 2 

2 .2 ft "H t 



/c«i 



2t 



The second limit (122) becomes 

1 _ e -I(3A A + KA )i 3A A + K A 



lim 
Vi 



-Ak-h 



3A A + K X 



Ax (3k h + \/l - \k h )(Sk h + 3\/l - Ak w ) 



2i 

3k« 
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and consequently (120) follows. To proceed with the proof of (114), we rearrange 

2 



,- dU x\ t ) 

lim £— — = Inn 

xyi oX \/»i 



Aa - K\ 



Ax 



3A A + k x 

A A 
Aa - k\ 

Aa - n\ 



Aa 



= lim 



Aa 

1 - e -^(3AA+^A)t 

3A A + k x 



.1 p— |(AA+«x)t 
2 



(-k^ + 2A a k w ) 



A A 



3A A + K A V2 



* p -i(3A x + n x )t ( %K H _ 



2A7 



K-H 



Ax 



2A A Kn 

e -^(A A +K A )t _ g-A x t 
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By means of (117) to (119), the limit of the expression after the squared brackets in (123) becomes 
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and the limit of the expression in (124) becomes with (120) 
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By putting (125), (126), (127) together with (120) we finally end up with 
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which finishes the proof of Lemma A. 6. 
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Proof of Theorem 6.6 Recall from (59) the approximative offspring-distribution-parameter /3; 
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by starting from Theorem 6.1(a). Correspondingly, 
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Additionally there holds 
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For k.a > 0, we apply the upper part of formula (72) as well as (128), (129) to derive 
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(129) 



For k_4 = (and thus k-h > 0, /?7 = 1, aT" = i]/a 2 ), we apply the lower part of formula (72) as well as 
(128), (129) to obtain 
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Let us now calculate the "converse" double limit 
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This will be achieved by evaluating for each t > the two limits 
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which will turn out to coincide; the involved lower and upper bound D\ t , D^ t defined by (66) and (66) 
'x,t 



satisfy limAyi D% t = limAyi D\ t = 1 as an easy consequence of the limits (cf. 64) 
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(131) 



as well as the formulae (117), (118) for the case ka = 0. Accordingly, we compute 
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(133) 



Niels B. Kammerer and Wolfgang Stummer/ Distances of branching processes 



12 



For the case ka > 0, one can combine this with (131) and (111) to end up with 
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(134) 



For the case ka = 0, we continue the calculation (132) by rearranging terms and by employing the formulae 
(112), (113), (117), (118) as well as the obvious relation ^ - A=fA. = J- to obtain 
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Let us now turn to the second limit (130) for which we compute analogously to (132) 
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(136) 



For the case ka > 0, one can combine this with (131), (133) and (111) to end up with 
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(137) 



For the case ka = 0, we continue the calculation of (136) by rearranging terms and by employing the formulae 
(114), (117), (118) as well as the obvious relation lmiAyi -r — A 7x737* \ = — to obtain 
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Since (134) coincides with (137) and (135) coincides with (138), we have finished the proof. □ 
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